something/someone which/who is the 'object' of 
measurement) become here crucially important. 

As a matter of fact any kind of measurement requires 
a standard experimental condition in which 
observations are made and the system's parameters are 
adjusted. This standard condition is a fixed frame of 
reference, an 'observer' who (from that classical 
perspective) is supposed to be 'impartial' and 'external 
to the system. 'Omnisciency'is often ascribed to him/her, 
so that all the outcomes that he/she regards as relevant 
to his/her enquiry are known beforehand. In these 
conditions the observer is also (supposed to be) able to 
entertain (deterministic and/or probabilistic) 
hypotheses about the occurence of such outcomes (one of 
which, stemming directly from the mechanicist 
influence upon the present-day vision, is based on the 
idea that whatever he/she is observing may be treated as 
a machine). As a further point, the 'system' which these 
outcomes are assumed to be 'expressions', i.e, the focus of 
the observer's attention (distinguished from the whole 
flux of events he/she is able to recognize due to some 
feature or property he/she deems relevant) is regarded 
as an 'impersonal entity' and unequivocally demarcated 
as an it, being in most cases divided into 'parts' (as 
'emmitter'"receiver', 'causative input'"output effect', 
etc.). One of these parts, the so-called 'environment', is 
designed by the experimenter/external observer to react 
in a particular way to the 'object/it' under scrutinity; 
failing this the experimenter acts as though he/she were 
the environment. In these conditions, however, his/her 
actions have to obey a strict and predetermined set of 
rules. On the other hand, and apart from these 
interactions (usually represented in terms of 
'independent' variables), all other fluctuactions (deemed 
irrelevant) are minimised; the idea here is to held the 
experimental conditions so that the experiment can be 
'exactly' replicated again and again. Experiments are, 
finally, caried out in such a way that (if possible) only 
discrete variations of one independent variable and the 
observation of one or few (observed or 'dependent”) 
variables are involved. 

Let me apply these tenets to the situation which Fig.l 
depicts. The picture comprises three main boxes dealing, 
respectively, with a real or experimental universe, an 
external observer and a modelling universe (for this 
observer). The experimental universe contains an 
'experiment' (meaning any constructible and observable 
situation) in which, for the sake of brevity, the observer 
has made a distinction between a 'subject'"organism' 
and a 'regulator' which, in this case, is a predefined (or 
programmed) 'environment'. This distintion 
organism/environment (the grounds of which are the 
same as those for making it in the modelling universe) is 
supposed to be held independently from the nature of 
subsequent observed interactions; the spatial 
boundaries there shown are, in a certain sense, a 
symbolic representations of this hypothesis. As a further 
assumption the observer is supposed to have a 
description of his/her model and a (possibly) imperfect 
description of the experiment. Obviously, the form and 
scope of these descriptions vary either from experiment 
to experiment or from observer to observer. But 
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whatever the experiment (or whoever the observer) may 
be, a mandatory requirement must be satisfied: Lhe 
model shall be communicable to ther observers using 
some kind of 'observational language' L. This language 
becomes a metalanguage 1.* when the observer talks 
about models as a class and also when he/she talks about 
real strutures and laboratory experiments. To avoid 
particular and detailed descriptions (meaningless in this 
context) let us represent descriptions in general by the 
symbol €, 

As previously emphasized the observer is endowed 
with the capability of operating upon the models and/or 
experiments under scrutinity. An operation that is 
unrestricted (apart from the rules of the observational 
language) is represented by a filled-in parametric arrow 
penetrating whatever is operating upon. "Tuning' or 
restricted operations (which symbolize adjusting 
parameters so that models fits the facts or adjusting an 
experimental desing so that it satisfies a model) are 
represented by a blank (parametric) arrow. 

Let us imagine that the experimenter's objective 
seeks the construction of a 'theory' of what is going on in 
the interaction organism/environment. At this point | 
shall remind you (according to the present-day 
interpretation of the classical paradigm) that a theory, 
within which an observer can make inferences, is a 
metaphor which (in the case of Fig.1) asserts a stricl 
analogy* between a relation coupling 'observer's model! 


standar | relation which, in Lhe simplest case, 
takes the form "À is to B as C is to D'. Let us assume for the sake of 
conciseness that À and B belongs to an universe of discourse, say 
U, and C and D to another universe, say Uy. Let F and G be also 
two relations carrying elements such as A (or C) in its domain 
into elements B (or D) in its codomain. Thus, an analogy exists 
insofar as there is a relation between F tin Uçjand G (in U,.. 
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for the organism' with a further relation coupling 'real 
organism' with 'special environment”. 

Tenure of this analogy is contingent upon its domain 
of interpretation which, ultimately, corresponds to an 
identification between the 'model universe' say U,, ard 
the 'experimental universe', say U,; it holds if the 
metalinguistic statement “there is an identification 
between U, and U,” is (regarded as) true. 

This identification (which, from the classical 
viewpoint, is usually obtained by specifying 
measurements and the like) is a stable condition in 
which the boundaries (or “boxes” in the particular case 
of Fig.1) that enclose 'model' and 'experiment' are 
supposed to retain their integrity and in which those in 
the 'model' are related to some or all of the 'experiment”. 
As a further, albeit mandatory, requirement variations 
which do not explicitly appear in the model are (assumed 
to be) 'constant conditions' in the experiment (the 
constants condition 'box' in Fig.1l has an unidirectional 
coupling directed to the organism). In other words, the 
material analogy underlying the theory has an universe 
in which these conditions are held constant as its 
domain of interpretation. 

In Fig.1 (as well as in pratically all the situations in 
which a physical experiment or an it is involved) the 
external observer carries out two kind of operations. On 
the one hand, by means of unrestricted operations 
he/she poses and tests hypotheses (typically, as 
previously stressed, a causal hypothesis about the 
organism includes some relation between “input 
stimuli” and “output responses”) which, if valid , allows 
him/her to make inferences within the teory. On the 
other, by means of 'tuning' operations, the observer can 
adjusts his/her model for the organism so that it can be 
identified under the measurement procedures he/she has 
at his/her disposal (which is tantamount to saying that 
tuning' secures the result of restricting the analogical 
universe in such a way that the analogy underlying the 
theory remains always strict). To this extent, contingent 
propositions may be (regarded as) true and the material 
analogy of the theory holds, in its proper domain of 
interpretation. 

These aspects being clarified let us assume that the 
experimenter/observer has in mind to bring up a 
'behaviour' (ultimately, a relation between the 
organism's 'responses' and the 'stimuli' to whish is being 
submitted). In order to achieve this goal (and according 
to the language of the classical paradigm) that part of 
the environment deemed relevant for the observer's 
purposes (the box labelled 'special environment to the 
hypotheses being tested") is suposed to act upon the 
organism in accordance with a predefined and rigid 
pattern. Since, however, this pattern is ultimately 
controlled by the experimenter, everything works as 
though the observer him/herself were causally acting 
upon the organism. This means that he/she will deliever 


" (conta) Ana Tógies are named “strict” if this relation is either 
an isomorphiam (one-to-one correspondence) or & homomorphism 
(many-to-one relation, preserving correspondence). But if F and 
G holds only in part of the analogice! universe (specified by an 
indefinite sequence of properties P, of U, and P, of U,) then the 
analogy is said to be locally strict at « given instant" (though it 
may become strict in further instanits) 


stimuli or 'reinforces' (perhaps contingent upon previous 
events) according to some hypotheses and for each 
stimulus the observed events (outputs, responses or 
'reactions') are described by the experimenter using the 
metalanguage L*. In order to enlarge his/her 
descriptions, the observer is still endowed (if desired and 
possible) with the capability of changing the 
environmental characteristis through the filled-in arrow 
previously referred to. Sometimes (as a further 
refinement) these transformations of the environments 
predefined pattern may be monitored by a sequential 
sampling system; but, in these conditions, the possible 
pattern must be known in advance (an hypothesis which 
stems directly from the assumptions on which the 
probabilistic calculus is based). Finally, the observer's 
experimental objectives are (supposed to be) achieved 
when a certain number (regarded as enough) of 
stimulus/responses or input/output data is reached, this 
in order to describe the organism (or, more correctly, a 
model of the organism) as a 'finite automaton' (in 
discrete mathematics) or by means of transfer functions 
or transfer matrices (always that such input/output 
relation are assumed to be continuously time-variable). 

Two aspects of this representation are crucial for the 
further understanding of the perspective from which the 
classical framework looks upon interactions in general: 
the observer's choice of a description and the limited 
scope of his/her parametric operations. Between them 
they determine the (model of the) environment as a 
finite-state automaton, a result which (together with 
what was asserted about the proper model of the 
organism as another automaton) lead to some further 
results. 

On the one hand, the interaction 
organism/environment becomes, thenceforth, 
interpreted in terms of an exchange of symbols in a 
formal language'. Despite its agreement with automata 
theory and mathematical psycholinguistics, this 
language' is, however, a syntatic entity only, never a 
language in its ordinary sense (in which, together with a 
syntax, a semantic and a pragmatics have also to be 
taken into account). Às a further consequence of this 
exclusion of semantics and pragmatics, the models of the 
organism and of the environment may be equivalently 
regarded as abstract automata (for accepting/generating 
strings of symbols) whose 'internal states' are the 'non- 
terminal" symbols of the formal language used and 
whose 'structures' can also be regarded as 'formal 
grammars', On the other hand, looking upon the system 
from the outside and depicting symbols or symbol-string 
occurences, the external observer becomes 
(theoretically, at least) able to construct contingency 
tables and/or to tabulate combinations of stimuli and 
responses. Next, by forming ratios of the frequencies 
according to which the combination of events occur, he 
may also estimate events probabilities, conjoint event 
probabilities and so forth. Briefly, he may build up the 
desired 'theory' of what is going on - this if the 'organism' 
under examination is, for example, a mouse, a pigeon or 
a reflex subsystem of a man. However, if the object of the 
observer's analysis is no longer an it but the human qua 
sentient being (say, a psychological entity where 
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everything emphasized about mentation has really to be 
considered) the construction of such a theory based upon 
the premises of the classical paradigm faces enormous 
and, even, unsurmountable difficulties. Because directly 
linked with the endeavoured extension of that paradigm 
to human psychology the following aspects are 
especially poignant. 

Firstly, it is clear that the meaning of such terms as 
'stimulus', 'parameters', 'responses', and so forth 
(derived from the language of physiology and 
behavioural psychology) becomes very strained when 
the mind's functioning is altogether considered. Not only 
because there are situations clearly out of court 
(starvation, for example, in terms of psychological work 
at least) but also because what really happens when this 
type of human interactions is involved (in which the 
existence of an experimenter or 'agent' and of an 
experiencing being or 'patient'is implicitly subsumed) is 
that the experimenter issues instructions in the belief 
that they will induce in the 'under experiment" human 
being an appropriate 'mental set". The trouble is that 
these instructions establish constancy (if they do so at 
all) not because they are imposed like a physiological 
constraint but - and this is the crucial aspect - because 
they are interpreted, understood and accepted. Because 
of this, instruction giving and instruction accepting are 
outside the scope of the classical paradigm; they belong 
to a normative framework rather than to functional 
situations bounded by conditions the experimenter is 
able to set up independently from the experiencing 
being. 

Although related to this criticism, the second aspect 
deepens (from a slightly different perspective) the 
automata theoretical approach which the fixed 
description that the classical paradigm underlies is the 
contemporary version. As | previously emphasized this 
fixed description determines the input and output states 
of the automata which organism and environment are 
supposed to be. But if this result is added to the 
experimenter's parametric restrictions another 
conclusion may be stressed out: that the environment's 
internal states cannot be changed as a consequence of 
any act (like attending to something new) that the 
'organism' might perform outside the limits of the initial 
description. In these conditions it is therefore possible to 
apply to the interactive scheme of Fig.1 the argument of 
Heinz Von Foester, who demonstrates that any 
behavioural sequence in that representation is, in some 
way, predetermined. In this sense the proposition "the 
coupled system organism/environment undergoes 
adaptation” becomes a pure tautology; the system must 
do so unless the entire paradigm is misapplied (for 
instance, because some condition is not held constant as 
required by that framework). In consequence, all modes 
of adaptation are predictable and none of them can 
justify (for example) not only attention, evolution or 
learning but also that crucial (and dynamic) feature of 
real human interactions of which the existence of 
'autonomous' (or 'operant') responses is the foremost 
externalized expression. For the classical observer these 
responses must, usually, be minimised (because 
irrelevant' for his/her confined purposes); however, 
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dealing with human qua sentient beings they, willy 
nilly, exist and the regulating/unregulating process they 
underly is a matter which, albeit hardly explicable from 
the pure classical approach, has factually to be taken 
into account. 

The third of the aforementioned difficulties comes up 
with this mention to operant responses. As a matter of 
fact, since such responses are commonly referred to other 
than 'non-external' sources of stimulation, then the 
(paradigmatic) distintion between 'objective' and 
'subjective' as well as of its cognate emphasis upon 
'objectivity' alone are both placed in a rather awkward 
position. Let me justify this assertion from a more 
detailed perspective. 

When applied to the classical experimenter/observer, 
'objectivity' usually underlies two further assumptions. 
On the one hand, he/she is (assumed to be) 'specialized' 
in some area, i.e he/she aims the discovery of the 
'invariant truth' about the activity or behaviour of a 
particular type of system; this in such a way that no 
personal idiosyncrasy is expected to influence his/her 
judgements. As such, this specialized observer must, on 
the other hand, use a specific 'reference frame', i.e, a 
body of knowledge in which statements are related to a 
common, albeit 'scientific', language; this reference 
frame determines the kind of enquiry that he/she 
regards as relevant and, consequently, the set of 
physical attributes (variables, parameters and so forth 
together with their numerical values) that he/she must 
select in order to achieve his/her previous objectives. 
However, even in the realm of physical systems, this 
apparently cogent picture conceals criticisms. 
Hypothetical complex systems are hardly related to 
their observed activities/behaviours and in simpler 
system the specialized observer is always limited by the 
number of states he/she can usefully distinguish in any 
experiment (which, putting the matter differently, is 
tantamount to saying that the experimenter is generally 
'uncertain' about the activities of his/her hypothesized 
system). One of processes used for decreasing this 
uncertainty, consists in the replication of experiments, 
each of which validates only a sub-system of the 
hypothetical system he/she assumes to exist. But, even 
in these conditions, the observer has to assimilate 
results, a fact which will imply, in turn, that the 
experimental observations must be not only comparable 
(remember what was asserted about the 'constant 
conditions' of organism and environment in the situation 
depicted in Fig.1) but (if this is the case) also 
transformed according to the composition rules which 
the (chosen) reference frame presupposes. 

Distinct interpretations of the (hypothetical) system 
under scrutinity in the light of its activities/behaviours 
are therefore possible whenever some of the following 
situations occur (others will not be referred to here): the 
chosen reference frame or its composition rules are 
changed for any reason whatever (history of sciences is, 
for instance, a good example of these changes), the 
attributes deemed relevant in different frames are either 
not-measurable in a physical sense (which, to a great 
extent, is the situation where thinking, as mind activity, 
is involved) or, if measurable, prove to be incomparable 
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(which implies that results from experiments in these 
frames are not assimilable and, consequently, that the 
language' used varies from observer to observer, raises 
severe ambiguities or, even, engenders different types of 
truths”. 

If these comments are added to what was asserted 
about the 'constancy of conditions' in organism and 
environment ('automous responses' included) then it is 
almost needless to remark how these situations beset 
not only the distinction between 'hard' and 'soft' sciences 
I have been mooting so far but also the question 
concerning real and free (syntatic/pragmatic) 
interactions amongst human beings - especially when 
one of them (or both) bear(s) in mind the interpretation 
(and further understanding) of the others ultimate 
thought processes. 


3 Summary 


Humanistic systems (treated by 'soft' disciplines such 
as sociology, psychology, education, anthropology, etc.) 
deal with a host of puzzling questions hard to overcome 
by means of the mental framework which, since the 50's 
in special, has ruled the broad mainstream of science. 
One of the major features of this way of thinking (which, 
for short, was named 'classical paradigm') as well as of 
the (cognate) factors underlying its limitations in the 
analysis of such humanistic system were briefly 
surveyed in the foregoing considerations. Amongst them 
especial emphasis was laid on the widely extent belief 
(consonant with that classical paradigm) which looks 
upon the human being not as qua sentient entity but 
rather as an it (in reality as a finite-state machine) 
whose behaviour (individually or socially taken) is 
(supposed to be) described in the syntatic (albeit 
semantically and pragmatically neutral) language 
either of 'stimuli'/'response' (physiology or 
pshychological behaviourism) or of 'input'/output' 
(General Systems Theory). From the major criticisms 
which this point of view begets, the emphasis laid on 
cognitive processes alone, with exclusion of (pratically) 
everything related to what human beings have of 
subjective, idiosyncratic and distinctive, the scanty 
identification of 'thinking' with (a rather restrictive 
notion of) 'reasoning' (based on the timeless and 
supposedely 'generaly shared laws' of classical logic), the 
strained meaning which terms such as those of 
input/stimulus' and 'output'/ 'response' (the former in 
particular) acquire when applied to the human mind's 
functioning (altogether considered), the exclusion 
(because deeimed irrelevant) of 'autonomous responses' 
and the features ascribed to the (classical) 'external 
observer' ('objectivity', 'omnisciency' and so forth) were 
(some) of the aspects more closely examined. Especial 
relevance was laid on all the situations in which any 
social interation is involved (even in the most 'exact', 
'scientific', discussions). This because each one of the 
interlocutors engaged in such a conversation have 
always to interpret, understand and (partially, at least) 
accept the others' utterances - a situation which does not 
pertain to the classical paradigm but, rather, to 
normative and relativistic frameworks. 


The analysys of these frameworks as well as of the 
matters they deal with will be develved in the following 
monograph. 
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RESUMO 


Neste artigo, através da crítica da visão 
paradigmática da termodinâmica, pretende-se 
evidenciar uma nova perspectiva que se adquiriu nesta 
matéria. A generalização da conceptualização cinética 
de calor (calor como energia interna) permite a 
generalização do postulado de Kelvin: "É impossível, 
sem outro efeito, converter energia interna em 
trabalho" (não se refere calor ou reservatório de calor). 


Introduction 


The conceptual foundations of thermodynamics 
present subtle and interesting difficulties. The 
conceptual works of Caratheodory [1] and Callendar [2] 
are good examples of that at the beginning of the XX 
century. A more recent work defines the generalized 
view that has survived to this day - the interesting and 
well known book of E. Fermi, Thermodynamies [3]. We 
call that view the paradigmatic view of 
thermodynamics. In this paper the method to introduce 
and explain the new approach proposed is related to the 
above mentioned book in so far as its fundamental 
concepts are criticized and the book can therefore be 
regarded as a paradigm, although with interesting 
peculiarities. (In the text we refer to Fermi's book as 
(F.p.x), giving the number of the page x). In this sense 
we are not going to criticize especially E. Fermi's book 
but the paradigmatic physical approach of the middle 
XX century. This has the obvious advantage of making 
it easier to pass from that view to the new approach 
proposed, with simultaneous corroboration. The use of a 
well known book with physical interpretations, like 
Fermi's, permits to achieve this goal. In fact, Fermi 
sacrifices the logic consistency of the phenomenological 
view to the obvious interpretation emerging from the 
kinetic view. 

In point 1 it is questioned if thermodynamics is an 
autonomous physical domain supported by the heat 
concept. This is a fundamental question because 
between two equilibrium points the existence or non 
existence of heat is dependent on the concept of heat in 
association with the concept of System [4]. In fact we can 
separate the energy interaction terms into a work term 
(gravitational potentional energy change or equivalent) 
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ABSTRACT 


In this paper, through a criticism of the 
paradigmatic view of thermodynamics, we uim at 
showing a new perspective attained in this matter. 
The generalization of heat as internal energy 
(generalization of the kinetic energy concept of heat) 
permits the generalization of the Kelvin postulate: "It 
is impossible, without another effect, to convert 
internal energy into work" (no reference to heat or to a 
heat reservoir). 


and an internal energy (Heat 1) variation term. The 
energy exchange between two parts of the system may be 
classified as heat exchange (Heat II). This heat 
(exchange) is not equivalent to the classification of the 
internal energy term as heat. We can have an exchange 
between work and internal energy and we can say that 
the internal energy of the system increases by a quantity 
equal to the work energy term - the work is transformed 
into heat. The heat exchange between two subsystems 
(Heat II) is not equal, of course, to the internal energy 
increase of each one of the subsystems because of the 
existence of the work term. 

If we use the word heat as internal energy (Heat I ) it 
is possible to generalize the ordinary enunciations of the 
"Second" Principle of Thermodynamics and also 
eliminate the essencial character apparently permitted 
by the "First" Principle of Thermodynamics. 

The use of the word heat has been recently (once 
again) proposed for the quantity entropy (Heat II) [5]. If 
we use the word heat for the internal energy we have 
most of the properties pointed out by H. Fuchs [5] with 
no need to affirm that we have a heat increase in a free 
expansion of a gas and the word heat can obviously be 
used in the energetic sense emerging from the 
fundamental work of Joule [6]. 

Terminological conflicts are also in the origin of the 
mistake pointed out in 2 and called the paradigmatic 
error of Thermodynamics. 

In 2.1 we begin by analising the confusion between 
the concept of quasi-static transformation and reversible 
transformation. 

In 2.2 we clarify the notion of work and, particulary, 
the notion of work for a reversible transformation. 

In 2.3 the problems analysed in 2.1 and 2.2 are related 
to the energy conservation principle and the concepts of 
heat - the meaning of the "First" Principle of 
Thermodynamics is analysed and three examples of the 
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conceptual difficulties emerging from the paradigmatic 
view are refered to at 2.3.1,2.3.2 and 2.3.3. 

In 2.3.1 the meanings of specific heats are related to 
Energy Conservation Principle avoiding the 
paradigmatic error. 

Since the concept of temperature can be derived from 
energy and entropy [7] and related to the ideal gas [8, 9], 
it is possible to consider temperature a derived and not 
essential concept and also eliminate the tautological 
association between the Kelvin temperature and the 
ideal gas absolute temperature. This is refered to in 
2.3.2 and developed in the article of reference [25]. 

In 2.3.3 we relate the notion of adiabatic 
transformation to the entropy variation and to the 
“First” Principle. 

In 2.3.4 we deal in detail with the meanings of heat 
and propose the use of the word heat as internal energy. 
For a transformation in contact with a heat reservoir the 
heat exchange (Heat Il) corresponds to the heat 
variation of the reservoir (Heat |). The essential and 
general character apparently given by the formal 
expression of the Ist Principle to the quantity dQ, 
different from the internal energy, is eliminated. In fact, 
the use of the formal expression of the first principle 
together with the confusion between the concept of 
quasi-static transformation and reversible 
transformation originates the paradigmatic error 
analysed at 2 and leads to another "subtle error" 
considered in another article and refered to in 2.3. 

In 3 a new perspective is acquired for the "Second" 
Principle from the internal energy conceptualization of 
heat. In fact the entropy concept can be derived from the 
equilibrium tendency [10] and from the asymmetric 
heat variation [11]. Instead of a Ist and a 2nd Principle 
we propose an energy-entropy principle [11]. 

In 4 the energy-entropy principle is interpreted 
microscopically in a straightforward and covenient way. 

In the text the words system and subsystem are used. 
The System is the ensemble of all subsystems that 
intervene at the transformation. This being so, the 
transformation is always adiabatic for the system. 
Therefore, for the system, it is tautological to say that 
the transformation is adiabatic |4]. 


1 Whatis Thermodynamics? 


In the introduction of his book Fermi tries to define 
the thermodynamics domain as the domain of heat. 
Curiously and paradoxically, this introduction displays 
the ordinary criticism to Carnot's conceptualization of 
heat. This generalized criticism is paradoxical because, 
as is well known and accepted, Carnot's work is the 
origin of the essential concept of thermodynamics - the 
entropy concept (see Appendix). 

As we will see, the phenomenological interpretation 
of heat is apparently inconsistent with other 
interpretations (the kinetic energy interpretation or the 
generalization of it - the internal energy interpretation 
proposed). 

Although Fermi considers thermodynamics (F.p.15) a 
heat science, in this introduction he says that 
thermodynamics may be considered a special branch of 
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mechanices [10], the statistical mechanics (F.p.IX). Then 
if we think that thermodynamics has a non mechanical 
concept, like heat, as Fermi claims, the paradox is 
evident. This false conflict between mechanics and 
thermodynamics, if accepted, as it is still accepted today 
together with other and convenient confusions, makes 
thermodynamics an apparently mysterious physical 
subject. The mystery however, can be avoided. 


2 The Paradigmatic Error 


2.1 The confusion between the concept 
of quasi-static transformation and 


reversible transformation 


One of the basic concepts of thermodynamies is the 
concept of reversible transformation as a transformation 
with intermediate equilibrium points (F.p.4). This 
conceptualization is the origin of what we call the 
paradigmatic error [12, 13]. One example is enough to 
characterize the problem. In fact if we idealize a free 
expansion of a gas with intermediate equilibrium points 
[14], it is impossible to make the gas return to its initial 
condition. This being so, the transformation satisfies the 
equilibrium condition but nevertheless it is impossible 
for the gas to return to the initial state - the 
transformation is irreversible. However, if we consider a 
system composed of the gas and a thermal reservoir, we 
can apparently have a return to the initial state (after 
free expansion the gas is compressed and forced back to 
the initial conditions) as long as the thermal reservoir is 
big enough to hide the impossibility for the actual 
system (subsystem gas + subsystem reservoir) to return 
to the initial state (only in the limit when the mass of the 
reservoir is infinite can the gas return to the initial state 
but even in this case the reservoir itself cannot) because 
the energy of the reservoir is altered. Although the gas, 
after free expansion, has returned to its initial state due 
to the compression in contact with the thermal reservoir, 
this reservoir retained the energy corresponding to the 
work done during the compression. The reservoir may 
obviously exchange energy with another subsystem and 
return to its initial state. But the new system in contact 
with the reservoir will, also obviously, be prevented from 
returning to its initial state, unless it enters into contact 
with another subsystem. And so on and on (see fig. 1). 

For the system (association of subsystems) the 
transformation is irreversible because the system will 
never be able to go back to its initial state. Because of 
what has been said, it is essential to distinguish between 
the whole (System) and the parts (subsystems). In the 
text the word system is used in this global sense, not in 
the common, currently used sense, of subsystem (the 
word System, with capital S, has been used to call 
attention to this). Only if the exterior force is equal to 
the interior force and does a work [15, 16] (with an 
energetic meaning) permitting the return to the initial 
conditions do we have, of course, a reversible 
transformation (see 2.2). 

To affirm Lhe quasi-static conditions is not enough to 
define the reversible condition. 
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2.2 Thework expression for a reversible 


transformation 


The energy principle and the tendency to equilibrium 
permits to write the expression (F.p.11) 


W=AU (1) 
where W is the work (associated with the rise or fall of a 
weight) [17, 18] and 4U is the System internal energy 
change between two equilibrium points (Fermi used for 
work the symbol L and the sign convention is such that 
L=-W). 

For a reversible transformation we can write 
dW=-pdV where dW is the infinitesimal work and p 
and dV are the equilibrium pressure and the 
infinitesimal volume change. In this case, and only in 
this case, can we write 

dW= —pdV = dU (2) 


Actually, the internal energy is also an entropy 
function (U =U(V, 8))[11]. Then 


aU FTA 
au=( 5) ave (S) ds. (3) 
av ), as ), 
So 
dU = -pdV+TdS (4) 


where T=(9U/98), is the Kelvin temperature and p and 
S are respectively the pressure and the entropy. 

For a reversible transformation the infinitesimal 
entropy change of the System is zero (dS = 0). 

Hwe have an infinitesimal irreversible change 


dU = -pdV+TdS (5) 
with dS >0, then 
dW=dU=-pdV+TdS 


and 


dW < —pdvV. 

For an irreversible infinitesimal change the work is 
not equal to -pdV. 

The paradigmatic error is to write dW =-p dV for all 
quasi-static transformations. Only for the reversible 
transformations this is correct. Formally, of course, we 
can write dW=-pdV but if we insist on calling this 
expression the work term we are comitting what we call 
the paradigmatic error [12, 13). 

2.3  The”"First” Principle of 
Thermodynamics 


Fermi (F.p.4) gives to the reversible transformation 
the quasi-static definition (see 2.1). Only for the 
reversible transformation, as we saw (2.2), is the work 
term expressed by -pdV. This is a source of confusion 
because we have an expression and an operational 
means to calculate W=[-pdV, a term which does not 
generally correspond to the energetic work. 

In chapter 2 (F.p.11), Fermi says that the First 
Principle is essentially the energy conservation 
principle for "thermodynamic systems”. If the 1st 
Principle is only the energy conservation principle there 
is no problem. But this is not so if the Ist Principle 
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introduces a new quantity dQ. In fact Fermi begins by 
introducing the internal energy concept for "mechanical 
systems”, as he says (F.p.13). For "mechanical systems” 
he concludes that W=4AU and he notes that the existence 
of energy imposes that W (the work term) be 
independent from the trajectory (F.p.12). But in two 
pages (F.p.14) he says that this is not so because he 
thinks that thermodynamics belongs to a 
non-mechanical domain characterized by the heat 
exchange dQ (F.p.14 and F.p.15). (The exact sense of the 
word system has been previously explained in order to 
avoid such classifications, 1.e., mechanical and 
thermodynamical. In this way, System corresponds to 
what Fermi designates as “mechanical system", and 
subsystem to what he designates as "thermodynamical 
system”). 

IH we have a System with a great number of particles 
we can have a tendency to equilibrium if we change the 
exterior conditions [10]. Then we can write between 
two equilibriums points 


W=AU=AU+AU, (7) 


for a System composed of the subsystems i (interior) plus 
a thermal reservoir (subsystem e (exterior)). 
Of course we can also write 


AU = W+Q (8) 
and have 
= AU. (9) 


Curiously enough Fermi (F.p.29) also uses the word 
heat in the sense of internal energy and this is the good 
terminology proposed. In fact if we aim at analysing 
with generality the means of work generalising the 
concepts acquired at simple "mechanical" situations, we 
have difficulties. For variable mass systems the force 
conceptualization as f=dp/dt and the formal use of the 
Ist Principle doesn't permit an universal interpretation 
of the quantities involved [19, 20]. For subsystems in 
interaction through a movable wall, the separation of 
the energy interaction and the use of the First Principle, 
have originated "subtle errors" [21-24]. We shall 
therefore use quantity Q only in the sense of energy 
exchanged with a heat reservoir (Q=-A4U,). If we want 
to continue with quantities with physical meaning we 
have to be careful (see 2.3.4). 


2.3.1 The meaning of the specific heats 
Fermi (F.p.20) writes for an infinitesimal 
transformation introducing the specific heats C, and Cy 


dU — dW = dQ 0) 


and 


dO +pdi = dó: (1) 

Eq. (11) is not correct for quasi-static irreversible 

transformations (dW=-pdV for reversible trans- 

formations only). We can change the state of a simple 
System and write for an irreversible transformation 


dU=-pdV+TdS = dW. 
IfdV=0 


dU=TdS=dW=C,dT (12) 


and therefore 
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ey-n(6), (E), E 
Port Ny dr 


We can say that Cy is the energy needed to obtain a 
unit temperature change for constant volume. 
If P,., (exterior pressure) is constant then we have 


dU=-pdV+TdS = dW =dW+dW' = 
=CdT+(-P dV) 
p ext 


For Pe=p (p is the constant pressure at any 
equilibrium point), 


(13) 


(14) 


CdT=TdS =dw. (15) 
Then 
Cc (5) (16) 
p lar) dr 


corresponding to the energy necessary to raise Lhe 
temperature by a unit at constant pressure. 

The experimental methods to determine €, and Cy 
are, in most cases, conceived by measuring a work term. 
Of course this energy can be measured with calorimetric 
methods and in this case corresponds to an exchange 
dQ=-dU,. But here we emphasize that the 
measurements of C, and Cy are made in physical 
conditions that are associated with energy measurement 
changes of the exterior subsystem (dU,) or a work term 
dW. Of course for the calorimetric methods these 
quantities can be associated with a heat exchange under 
such conditions the heat exchange has a clear physical 
meaning... 


2.3.2 The concept of temperature derived 
from energy and entropy and the 
ideal gas 

If we introduce the temperature concept from energy 
and entropy we obtain T=(dU/9S).. 

For an “ideal gas” Fermi introduces the energy 
dependence exclusively on temperature from an 
inexistent experiment (curiously enough Fermi says 
(F.p.22) that the temperature change is small for the 
Joule experiment ! ) and uses this condition to identify 
the Kelvin temperature (T=9U/9S) with the ideal gas 
absolute temperature (pV=N K A; A is the absolute 
temperature). This is tautological, as Fermi says 
(F.p.62), and can be avoided, but here only the problem 
is referred to [25, 26]. 

2.3.5 Adiabatic transformations 

On page 25 Fermi affirms that an adiabatic 
transformation is a reversible transformation for a 
thermally isolated system (dQ =0). But if we are not 


careful we can commit the paradigmatic error. In fact, if 
we write 


dU = dW + dQ 
with dQ =0 (adiabatic condition - without heat exchange 
(Heat I1)), and if we also write 
dU=-pdV+TdsS 
we commit the paradigmatic error if we identify 


and 


36 


dQ = TdsS. (8) 

The adiabatic condition dQ =0 only imposes dS =0 for 
a reversible transformation (note that we are not saying 
that Fermi commits the error of affirming that an 
adiabatic transformation is equivalent to a reversible 
adiabatic transformation, because here, on page 25, 
Fermi clearly says that he uses the word adiabatic as 
reversible adiabatic). 

Fermi (F.p.26) uses the equation pV* imposing the 
adiabatic “reversible” (quasi-static) condition on the 
expansion of an atmospheric mass gas. He affirms the 
poor thermal conductivity for the air but says nothing 
about the reversible condition. He commits the 
paradigmatic error as he confuses reversible with 
quasi-static. 

An isentropic transformation (dS =0) for a gas with 
state equation p=a u (a is a constant and u is the energy 
density) satisfies the equation pVa +! constant. For a non 
relativistic monoatomic gas a=2/3 (for a photon gas 
a= 1/3). For normal (p, T) condition this gas approaches 
the classical Maxwell-Boltzman behaviour and satisfies 
an equation pV=BT with B constant. We have a 
classical ideal gas (the photon gas exists together with 
the classical gas, also satisfying an equation pV=BT but 
B is not constant) [25, 26]. 

For an ideal gas (p=a u) [8,25], dS=0 implies 
pVa+! =constant [25]. For an adiabatic irreversible 
transformation at constant exterior pressure (equal to 
the interior final pressure) between two equilibrium 
points, the pressure and volume for the initial and final 
points can be fitted by an equation pVê=constant with 
óza+l, 

If the final pressure p, is not too different from the 
initial pressure p (p, =p,), thenó=a+1. 

Nevertheless the paradigmatic error is committed for 
the reason that the equation for the irreversible 
adiabatic equation is pVô for the two points and not 
pVa+! =constant for all the points (dS=0). Although for 
the reason previously stated, 1.e., the proximity of points 
l and 2, ó value is approximately a + 1. 


2.3.4. Two meanings of Heat 

On pg. 29,31 and 56 Fermi astonishingly uses the 
word heat in a different sense from that he carefully 
defines on pg. 17 without explicitly making the 
distinction between these senses. He continues to use the 
word heat only. This is another source of confusion. In 
fact, on pg. 56 Fermi explicitly affirms the production of 
heat by friction ("since the heat comes from work and not 
from another part of the system..."). The origin of this 
mistake is the fluidlike conceptualization of heat (if 
water (H O) passes from one reservoir to another, the 
chemical properties in the passage are still those 
associated with HO and the mass of H,O is also 


conserved but if, for instance, the water runs a mill, the 


water energy is not conserved, Even for a fluid there are 
properties that cannot be conserved - see Appendix). The 
"heat exchange” (Heat II) between two subsystems is not 
equivalent to the subsystems heat variation - the 
energetic conceptualization of heat doesn't permit a 
fluidlike character for heat. In this sense we must avoid 
the terminology heat exchange (in the sense of 
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indestructible substance with intrinsic attributes as 
H,O has) but we can use the expression in the sense of 
energy exchange between two subsystems in special 
conditions. 

This being so, we propose the following nomenclature. 
The word heat is used as synonymous with internal 
energy. We can only have equality for the energy 
exchange between two subsystems and the internal 
energy change of the subsystems if there is no work. But 
as we had seen, when there exists work we can identify 
dQ with dU . Then, also in that situation (contact with a 
heat reservoir), we can think in terms of internal energy 
change of the exterior subsystem (the heat reservoir). 
This being so, we can say that we have “heat exchange” 
(Heat II) when we are describing an energetic exchange 
between subsystems and one of these subsystems is a 
"source of heat", and heat variation (Heat | variation) 
when we refer to the internal energy variation. The heat 
variation of the heat reservoir is equal to the heat 
exchange (with minus sign), the Heat | variation of the 
heat reservoir is equal to Heat II! 

À rigorous and general description of the dynamical 
energy change is attained by the energy impulsion 
conservation principle but in special situations we can 
use with physical meaning the expressions 
dW=dU+dU =dU-dQ and use the word heat in the 
sense I and IÍ (in the sense of internal energy of the 
“system” and in the sense of internal energy variation of 
the heat reservoir). This clear physical interpretation 
led us to the energy-entropy principle considered below. 


3 The "fractal-like” character of the 
energy-entropy principle - the 
generalization of Kelvin postulate 


Fermi (F.p.30) gives one of the traditional 
enunciations of the “second” principle of 
thermodynamics the Kelvin enunciation: 

"It is impossible, without another effect, to convert 
heat, extracted from a heat source into work”. 

Fermi (F.p.29) says that work can be transformed into 
heat but the reciprocal is not true. When he explains 
that “a body can always be heated by friction receiving a 
quantity of heat (Heat | variation) equalizing the work 
done, he is clearly using the word heat with a different 
meaning from that he had previously and carefully 
constructed (F.p.17) and (F.p.56). He uses the word heat 
with meaning I, in the sense of internal energy. 

We can generalize the Kelvin postulate affirming 
that it is impossible to transform internal energy into 
work without other effects (no reference to heat or a heat 
reservoir). The Kelvin postulate is a particular case of 
this because it affirms that it is impossible to convert the 
internal energy of a heat reservoir into work without 
another effect. 

As an exemple, consider a pas in a cylinder with a 
movable piston. It is possible to transform the heat of the 
gas (Heat 1) into work by increasing the volume. But if 
the volume returns to the initial value, then the internal 
energy is bigger or equal to the initial value (the equal 
value corresponds, of course, to the reversible 
transformation). For the same volume there is work 
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transformed into heat. As a consequence of this the 
internal energy cannot be a function of volume. 
Therefore we can introduce the entropy variable (see eq. 
(19)) (Note that p=-(9U/9V)s and T=(9U/9S)y are 
derived from U=U(V,S) - through the Legendre 
transformations) [7, 11,25, 26]. 

It is intersting to note that if we have an hypothetical 
natural source of "movable cylinder-piston-gas 
apparatus”" (the evolution after the big-bang has 
originated such things...) then we can lift weights (we 
assume, of course, that the natural initial condition of 
the gas inside the cylinder has a pressure bigger than 
the atmospheric pressure) and with these weights at a 
higher level we can construct a hot source which is 
continually renewed by the "movable cylinder-piston-- 
gas apparatus” (this is analogous to the burn of a fossil 
combustible). Or with this apparatus it is possible to 
obtain electrical energy (for example) without 
constructing a hot source, as we do in an internal 
combustion engine. In this sense it is not true that two 
heat reservoirs are necessary to produce work[27, 28] 
because initially we have only one temperature. 

For a system we must conclude the following: 

The transformation of internal energy into work is 
impossible without other alterations (obviously different 
from internal energy variation, or derived quantities!). 
We can say that only transformations with increasing 
entropy are possible (in fact this affirmation is derived 
from the previous enuntiation). But as we had seen, we 
have generalized it we can transform heat for all the 
processes corresponding to an increase of entropy (of 
course, for a heat reservoir it is impossible to transform 
the reservoir internal energy into work without other 
effects and this being so, beginning with only one 
temperature, a transformation with increasing entropy 
is necessary to have a global positive work). It should be 
noted that the previous statements cause energy to be a 
function of volume and of entropy [11]. And since there 
can only be transformation of work into internal energy 
for the same volume, entropy variation should have only 
one sign (the positive sign was chosen). Therefore it is 
clear that a Carnot Cycle with global zero entropy 
variation does not contradict the previous affirmation 
because we must have an irreversible process to “"renew”" 
the heat reservoir. But in a more restrictive way we can 
say that only processes without negative entropy 
variation are possible. 

This corresponds to what we propose to call the 
energy-entropy principle: 

The internal energy is a function of entropy satisfying 
the following relation: 


UV, S,)= U(V,S,). 


The reversible transformation corresponds to dS = 0. 
We can chose (9U/98)>0. Then, dS>0. Ifso, S, >8.. 
In fact for a System with 


(19) 


| | oU 
U=U(V,S), au=( >) av+( =) dS=dW. 
oV S as v 
For a reversible transformation 
dW = —pdvV. 


Then 
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dV = (— Jar) dS 
— [) = 

aV s as . V 
and 


Then 
Since dU =0 for V equal to the initial value 


(r=( ou ) 
Ss /y 
andifT>0, then dS >0. 

The entropy variation has two Lerms. À term 
associated with the volume, P/T dV, and the term dU/T. 
When dV = 0, dU >0 and dU <0 only if dV<0. (The 
microscopic interpretation considers these two terms in 
phase space (see point 4)). 

Finally we note the "fractal-like” character of the 
energy-entropy principle. To obtain work we must 
transform internal energy in order to return locally to 
the initial conditions. Then because this is impossible 
only for global entropy increasing processes, we 
necessarily have a subsystem having an entropy 
increase and for a finite system it is impossible to have 
the "machine" (a subsystem) return to its initial 
condition. If locally (where the machine is working) we 
have entropy levels no longer allowing the existence of 
the process, then we have to construct a new machine (a 
new process) permitting the decrease of the local 
entropy. 


dU=TdS 20 


lN mi nem Mi |] Ml 
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A triangle represents a cyclic process doing work W. The 
process can only be cyclic if there exists contact with the 
exterior represented by the rectangles and by the arrows. 
This infinite recurrence process is symbolically represented 
by sucessive triangles and rectangles. 


Figure 1 


But only a global increase of entropy is possible. Then 
we have the problem once again, corresponding to an 
infinite recurrence process, represented in fractal-like 
figure 1. With this global perspective we reveal the 
"fractal-like" character of the energy-entropy principle. 
This designation suggested by the fractal figures, points 
to an interpretation different from the "substance-like" 
interpretation arising from the paradigmaLic view 


4 The microscopic interpretation of 
the "second" principle 


Fermi, on pg. 56, and curiously enough, since he 
constructs the entropy change associated with the 
expression dS=dQ/T (only true for a reversible 
transformation, but as “reversible" is confused with 
"quasi-static" this expression is erroneously applied also 
to a quasi-static transformation), says that when a body 
is heated by friction it receives a positive quantity of 
heat and, because heat results from work and not from 
another part of the "system" (subsystem), the increase of 
entropy is not compensated by the decrease of entropy of 
the other part. The explanation is, of course, this one if 
we use heat in the sense of internal energy. But Fermi 
has defined heat with the meaning emerging from the 
First Principle. Without another explanation this is, of 
course, a source of confusion, 

The explanation is this: 

The work transformed into internal energy (Heat 1) 

dW = dU. 

For an infinitesimal transformation we write 

dW = dU = -pdV+TdsS. 


dV=o0, 
dW=dU=TAS. (20) 
Then 
W 
ms (21) 
TO T 


The expression is dS= dU!T (not dS= dQ/T ) and dU 
is the heat variation due to the friction work, indeed! 

All of this can be easily interpreted microscopically 
but the paradigmatic error (dS = dQ/T for a quasi-static 
irreversible transformation) with the confusion between 
Heat | and Heat II, originates interesting difficulties 
associated with the microscopic interpretation of the 
second principle. 

One of these difficulties is linked to the affirmation 
that the entropy of an "isolated system" can only 
increase (F.p.55). An "isolated system” is a System with 
constant energy (F.p.11). Here we have again a 
terminological problem if "isolated system" means 
“Lhermally isolated system” (in addition to this, there is 
the confusion associated with the concept of thermal - 
Heat | and Heat I1). Note that we are not saying that 
“isolated system" is equivalent to "thermally isolated 
system”, of course! (an “isolated system" is thermally 
isolated but the reciprocal cannot be true.) 

For a non isolated System as well as for an isolated 
System the entropy change is positive or zero. 
Nevertheless for a subsystem the entropy change can be 
negative. 

The entropy variation is due to the variation of the 
number of microstates. 

Only a transformation with an increasing number of 
microstates 1s possible. 

The energy contributes to this number of microstates 
and we can have a transformation of heat into work 
(with a decrease of Lhe associated number of microstates, 
the momentum space microstates) if there is an increase 
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of the total number of microstates. This corresponds to 
the energy- entropy principle. 


Conclusions 


In this paper we have shown that it is possible to 
eliminate the therminological conflits still existing in 
thermodynamis. 

Essentially, these conflits are the following: 

1. The word heat is usually associated with the 
expression of the First Principle dU = dW + dQ (U is 
the internal energy and dW and dQ are the infinitesimal 
work and the infinitesimal heat). The quantity dQ can 
be associated with the internal energy variation of a 
heat reservoir when the subsystem considered has an 
isothermal transformation (temperature constant). 
When the transformation is adiabatic (the energy 
exchange with the heat reservoir is zero), the internal 
energy of the subsystem considered changes due to the 
work energy term. For the system (association of the 
subsystems “system” + heat reservoir) the 
transformation is always adiabatic. Therefore the 
essential quantity of thermodynamics is the internal 
energy and not the heat exchange. It is possible to call 
the internal energy heat. But the internal energy 
(energy associated with the microscopic entities) is also 
an entropy function, the variable associated with the 
equilibrium tendency. The increase of entropy results 
from the increase of heat for the same deformation 
variable (it is only as a limit that energy and entropy do 
not change - this is the reversible transformation). Then 
instead of a Ist and 2nd Principle, Thermodynamics is 
characterized by a single principle, the energy-entropy 
principle (it is impossible for a machine to do work 
returning to the initial state for two non separate 
reasons: if the machine returns to the initial state the 
internal energy is the same and the entropy is also the 
same). But if the deformation variable is the same, the 
energy would have to be smaller for the work to be 
negative, which would correspond to a negative entropy 
variation. Only transformations with a positive entropy 
variation are possible. 

The word heat cannot be used simultaneously for the 
internal energy of a subsystem ("system") and for the 
energy exchanged between this subsystem and the heat 
reservoir because this energy is not equal to the internal 
energy variation of the “system” but rather equal to the 
internal variation of the heat reservoir. This is one of the 
terminological conflicts that has been solved [29]. 

2. The second conflict, although different, is 
intimately connected with the first. 

A quasi-static transformation is not equivalent to a 
reversible transformation. This confusion is related to 
the former in the following way: 

The "Ist" Principle 

dU = dW+dQ 
is erroneously related to the expression 
dU = —-pdV+TdsS 
derived from U = U (V,S). For an irreversible quasi- 
static transformation these two expressions can be 
formally related in the following way: 
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dW = —pdV 
and 

dQ = Tds. 
For a reversible transformation dS = 0, and then, dQ = 
O. For an irreversible quasistatic transformation 
dU + pdV 

T 

IH dV = 0 then dS = dUIT and not dS = dQ/T.k is 
important to note however, that, the paradigmatic error 
(ironically!) leads to the following: people think they are 
applying dS = dQ/T = dU /T + pdViT to a "reversible" 
transformation, when they are, in fact, applying it to an 
irreversible, quasi-static transformation (because the 
quasi-static transformation is erroneously called 
reversible). In this irreversible quasi-static 
transformation 


dU + pdV dQ 


ds = » but dS*& — 1 
já T 


Thus is understood and eliminated the conflict between 
the internal energy conceptualization of heat and the 
heat (exchanged) emerging from the “Ist" Principle. In 
fact the energy conservation principle preceeds the 
Clausius expression dU = dW + dQ. For a "mechanical 
system” (F.p.11) 
dU = dW. 

But we can write for a "system" (a " “ermodynamic 
system") in contact with a heat reservoir (this heat 
reservoir is a subsystem exterior to the "system”) 


dW = dU + dU.. 


Therefore 
dQ = —dU, 


and the "mechanical" condition is satisfied. 
The formal expression of the Ist Principle can induce 
thinking that the quantity dQ is the fundamental 
element of Thermodynamics. This is not so and 
originates subtle mistakes. In this sense it is important 
to eliminate the "Ist Principle”. This can be done with 
the construction of a single energy-entropy principle. 
The microscopic interpretation of the energy-entropy 
principle is straightforward and consistently related to 
the "phenomenological" approach proposed. The 
tendency to a new equilibrium point with bigger entropy 
results from a bigger volume or a bigger energy 
dS= — dU+ — dv. 
: à li 


The thermodynamic probability increases with the 
increase of volume or with the increase of energy we can 
say with the increase of the phase space volume. 


Appendix 


Carnot [30], in this book reasons by analogy with a 
machine that works at the expense of water passing from 
a higher into a lower level reservoir. 

If the water could, by itself, pass into a higher level 
reservoir, we would have a perpetual machine. Carnot 
associates the heat (or caloric) passing between two 
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sources at different temperatures with the water falling 
between the two levels. The quantity of water falling 
into the lower level reservoir is equal to the one coming 
out of the higher level reservoir. However, the energy of 
the water is not the same if the water has done work. 
There is, therefore, no energy conservation for the 
water. In the same way, there would be no energy 
conservation in the caloric - if the caloric existed! In fact 
the relevant property in this analysis is energy. So the 
analysis would be correct even if there existed a 
substance associated with the process! 

If the water could, by itself, (at the expense of its 
internal energy) pass from the lower into the higher 
level reservoir, there would exist a perpectual machine! 
Curiously enough, Carnot's analysis contains the 
energy-entropy principle! 
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RESUMO 


Este artigo descreve uma simulação computacional 
da Assinatura de um receptor digital na presença de 
desvanecimento selectivo. Descreve-se o método 
estatístico utilizado e apresentam-se Assinaturas de 
receptores não igualados e com igualação linear 
adaptativa, comparando as Assinaturas obtidas com 
uma aproximação analítica conhecida. 


1 Introdução 


O método de determinação da Assinatura de um 
receptor num sistema de feixes hertzianos digital 
apresentado em [1] utiliza um majorante da 
interferência intersimbólica (i. i. s). Admite que todos os 
símbolos interferentes têm amplitude máxima e que as 
interferências se somam em módulo. Seguidamente, 
determina por excesso, o número de símbolos 
interferentes, Ny, que leva o sistema a uma dada 
probabilidade de erro. Esta aproximação tem o mérito de 
tornar o método de cálculo da Assinatura independente 
das sequências de bits e permite obter uma forma 
analítica de resolver o problema. Parecem, no entanto, 
pertinentes as seguintes questões: Será o majorante 
suficientemente próximo da Assinatura real? Qual é a 
diferença entre a Assinatura real e o majorante obtido 
por via analítica? Para responder a estas questões, sem 
recorrer a métodos experimentais, ter-se-á de recorrer a 
simulação em computador. 

A simulação constitui uma forma de avaliar o 
desempenho de um sistema que não possa ser expresso 
numa forma analítica devido à sua complexidade e/ou ao 
envolvimento de variáveis aleatórias. Simula-se o 
modelo numérico do sistema durante um periodo 
suficientemente longo e estima-se as suas características 
com um dado intervalo de confiança. 

Neste caso concreto de avaliação de Assinaturas, à 
simulação tem a vantagem de permitir um melhor 
controlo das condições experimentais, uma vez que num 
sistema real o desvanecimento selectivo só se manifesta 
em períodos de tempo restritos. No entanto, a simulação 
requer recursos computacionais elevados para obter 
uma Ássinatura e, uma vez que estão envolvidos 
processos estocásticos no sistema, são necessárias várias 
amostras independentes do processo para conseguir 
obter as suas características como uma boa 
aproximação. Para conseguir a independência das 
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The present paper describes a computer simulation 
of a digital receiver's Signature, facing a selective 
fading. The statistical method used in the simulation 
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known analytic approximation. 


amostras e um intervalo de confiança suficientemante 
pequeno, para a precisão requerida, recorreu-se ao 
chamado procedimento por lotes (em Inglês 
"Batch-Mean". [2] 

Uma simulação com sucesso envolve algo mais do que 
construir um diagrama de fluxo do sistema em estudo, 
traduzindo-o em linguagem computacional. À teoria das 
probabilidades e das estatísticas são parte integrante do 
estudo das técnicas de simulação para escolher as 
distribuições de probabilidade das variáveis aleatórias 
geradas, validar o modelo de simulação e analisar 
estatísticamente os dados obtidos. 


2 Geração de variáveis aleatórias 


No sistema que pretendemos simular existe uma 
variável aleatória de entrada: a informação sob forma de 
bits. Estes, ao serem modulados em M-QAM! e enviados 
através de um canal de transmissão na presença de 
desvanecimento selectivo, chegam ao receptor e são 
recuperados com uma dada probabilidade de erro de bit, 
P.,, variável aleatória cujo valor médio depende dos 
parâmetros do canal b, fo e 1 definidos por Rummiler [3,4]. 
A Assinatura é obtida fixando 1, e para cada valor de fo, 
determinando o valor de b que leva a uma dada 
probabilidade de erro de bit, em geral no intervalo 
[10-6, 10-3]. Uma vez que o processo utiliza uma variável 
aleatória de entrada, os dados obtidos são também 
aleatórios. Para representar a informação a transmitir 
ter-se-á que gerar uma cadeia de bits equiprováveis e 
estatisticamente independentes. 

Um "bom" gerador de números aleatórios deve possuir 
algumas propriedades importantes: os números 
produzidos devem ter uma distribuição uniforme num 
dado intervalo (usualmente [0,1]) e não devem 
apresentar correlação apreciável. A geração de números 
aleatórios tem uma metodologia bastante complexa que 
se afasta do âmbito deste trabalho, recomendando-se 
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Averille Keltor [2] para uma apresentação completa dos 
vários métodos conhecidos de geração de variáveis 
aleatórias. 

Um dos tipos de geradores de variáveis aleatórias 
mais recentes são os chamados PMMLGC's (do Inglês 
“Prime Modulus Multiplicative Linear Congruential 
Generators"). Este tipo de geradores são vocacionados 
para uma execução rápida (em computador) e possuem 
um "bom" comportamento estatístico. A rotina utilizada 
vem descrita em [2] e funciona em máquinas baseadas 
em palavras de, pelo menos, 32 bits (incluindo o bit de 
sinal), assumindo que todos os inteiros desde -231 +] até 
231.1 são correctamente representados e calculados. Us 
parâmetros escolhidos permitem um comportamento 
estatístico óptimo e estão devidamente justificados em 
[2]. 


3 Intervalo de confiança 


A simulação do sistema não é executada por tempo 
indeterminado. Durante a simulação calcula-se um 
intervalo de confiança com uma precisão considerada 
aceitável e, quando esta for atingida, a simulação 
termina com uma estimativa das variáveis de saída, 
neste caso únicamente a probabilidade de erro de bit. 

Sejam Xy, Xo, ..., X, variáveis aleatórias 
uniformemente distribuidas, com média finita u e 
variância finita 02. Pretende-se construir um intervalo 
de confiança para yu e Lestar a hipótese complementar: 
H= Ho. 

Seja Z" a variável aleatória: 


[X (mn) (1) 
Vo? m 


e seja F,"(z) a função de distribuição de Z'“ para uma 
amostra de comprimento n, ou seja, F,"(z)=P/Z"' <=Z). 
O Teorema do Limite Central estabelece que: 


An) — 


FP (2) >q (2) 
onde p(z) é a função de distribuição de uma variável 
aleatória normalizada (u=0,02= 1) dada pela expressão: 


quando n> 


2) 


1 e = Ze 
& (z)= — | e “dy para —0< 2 <o ( 
V2n1 a 
Para n finito, o teorema mantém a validade 
substituindo o? por S2(n) 2 e a função de distribuição 
resultante é uma distribuição t-Student dada por: 
(nm) rw ara (3) 
=[X(n)-u]/V Sn)/n 
com n-1 graus de liberdade (n=2) e o intervalo de 
confiança a 100(1 -a)% para g é dado por: 
— o (4) 
AM Lt iai V S“(n)/n 


onde tn-1,1-a/2 é o denominado ponto crítico da 
distribuição t-Student com n-1 graus de liberdade. A 
estatística t-Student apresenta-se tabelada em [21]. 

A escolha de n deve permitir obter um intervalo de 
confiança com significância, ou seja, com correlação 
estatísticamente nula. 


E Estimador de Variância 
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A partir de uma amostra pretende-se testar a hipótese 
nula Ho: u=uy, onde 4, é um dado valor por hipótese. 
Escolhendo um valor para o intervalo de confiança 
(tipicamente 90-95%), faz-se o seguinte Leste: 


se | do | id aceita-se a hipótese H, 


n-1,1-a/2 


se | ár > rejeita-se a hipótese nula, 


n=1,1-0/2 
sendo a decisão tomada com 100(1 -a)% de confiança. 


4 Procedimento por lotes 


O procedimento por lotes começa por agrupar em lotes 
as amostras de uma variável aleatória, provenientes de 
um processo estocástico. À primeira vista, agrupar 
amostras não parece ter qualquer vantagem do ponto de 
vista estatístico. No entanto, demonstra-se [2], que tal 
procedimento reduz substancialmente o número de 
amostras necessárias para obter resultados com 
significado estatístico. 

Sejam Y, amostras de uma variável aleatória 
aproximadamente independentes com média E(Y)=v. 
Efectuando uma simulação com m concretizações e 
dividindo as várias observações de Y,,Y5,..., Ym em n 
lotes de comprimento [, assumindo m=n.[, então, o lote 1 
consiste nas observações Y;, Y9,..., Y,, 0 lote 2 nas 
observações Yi+1, Yiy2,-.Y21, etc. Seja Yjl), para 


j=12,...,n, a média das [ amostras no j-ésimo lote e seja: 
= Ea m 
F(n,D= > Y ADIn = > Y/m (5) 
É is | i=1 


a média da amostra total. Então, utilizando Y(n, [) como 
estimador de v, o intervalo de confiança constroi-se na 
forma: 

2 | 
14-88 Vit bull (6) 


Y(nD+t! 


onde: 


> Y, TD Y(n,1) 


Sy (n)=D——— 
é n—l 


(7) 


um 


sendo a expressão (6) análoga à expressão (4). 

Este método permite, com uma escolha apropriada de 
nel, obter um intervalo de confiança da variável 
aleatória a estimar. No entanto, persiste a escolha de n € 
!. Uma vez que a estatística t-Student conduz a 
resultados falseados se as amostras forem 
correlacionadas, é necessário escolher n e | por forma a 
que a correlação entre as amostras seja reduzida. Para 
esse efeito é necessário obter um estimador da 
correlação. O estimador apropriado é obtido através da 
expressão: |2] 

n-—l a 
à Ya tm Y 
% 
j= 


()— Y(n,0 (8) 


Dn d)= 
Ep (D)— Fa [) 


ia 
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Se o estimador indicar correlação (típicamente 
considera-se amostras incorrelacionadas se |p(n, )|< 0.3) 
ou se o intervalo de confiança não for suficientemente 
pequeno, então o conjunto de amostras não é 
significativo. Neste caso, calcula-se a média dos n lotes 
dois a dois obtendo-se n/2 novos lotes. Seguidamente, 
gera-se os restantes n/2 lotes, cada um com 2.l amostras, 
obtendo-se novamente n lotes à custa de mais m 
amostras. Repetindo os testes estatísticos de correlação e 
do intervalo de confiança verifica-se os novos resultados 
estatísticos, tornando o processo iterativo até que se 
verifiquem as condições estatísticas impostas. 


5 Simulação de um sistema M-QAM 


O sistema a simular é constituido pelo modulador 
M-QAM, o canal de transmissão na presença de 
desvanecimento selectivo e o receptor de M-QAM, nas 
condições apresentadas em [1]. 

Inicialmente é gerada uma sequência aleatória de bits 
(Cf. $2) que são convertidos em símbolos de M-QAM em 
código de Gray, obtendo-se assim os vários dibits ap; € 
Gqi- 

Seguidamente, é contabilizada a interferência 
intersimbólica dos vários símbolos nos instantes 
tn=to+n.T, (n=0, neZ), onde ty representa o instante 
ideal de amostragem defini... «em [1] e T, o período de 
símbolo em M-QAM. Obtém-se assim, as componentes 
em fase e quadratura dos símbolos na recepção, 
respectivamente, up(tg) e ug(to). 

Finalmente, no receptor, é tomada uma decisão acerca 
dos símbolos recebidos construindo-se uma sequência de 
bits na saída. Comparando-se esta sequência com a 
sequência gerada, avalia-se a probabilidade de erro de 
bit para um dado conjunto de parâmetros do canal. 
Pretende-se, fixando 1, e para cada valor de fy, obter o 
valor de b que leva o receptor a uma dada probabilidade 
de erro de bit no intervalo [10-8, 10-3]. 


5.1 Codificação de Gray em M-QAM 
Uma das vantagens da modulação M-QAM, dentro da 


classe geral dos códigos MAPSK 3, é apresentar um 
desempenho superior sem que a implementação seja 
demasiado complexa. Em relação à modulação N-PSK 4, 
também muito vulgar em sistemas de feixes hertzianos 
(note-se que 4-QAM é equivalente a 4-PSK), a 
modulação M-QAM é vantajosa para índices de 
modulação superiores a 8 devido à menor potência de 
transmissão requerida para assegurar a mesma 
probabilidade de erro [5). 

A modulação M-QAM, utilizando codificação de Gray, 
tem a vantagem adicional, em relação a uma codificação 
binária simples, de proporcionar uma melhoria do 
desempenho do sistema. Esta melhoria é obtida 
codificando os simbolos em M-QAM de forma a que os 
símbolos adjacentes só difiram de um único bit. Assim, 
existindo um erro na decisão, quer devido a ruído 
térmico, quer devido a i.i.s. ou à ambas as causas (ou 
outras), a probabilidade do símbolo resultante do erro 
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ser adjacente ao símbolo correcto é muito maior que a 
probabilidade complementar. Assim, no caso do símbolo 
errado ser adjacente ao símbolo correcto, utilizando 
codificação de Gray, erra-se apenas um único bit e em 
codificação binária poder-se-ão errar dois ou mais bits, 
dependendo do índice de modulação. 

A constelação de M-QAM em código de Gray não é 
única, existindo um conjunto limitado de possibilidades 
diferentes com o mesmo objectivo. No entanto, as várias 
possibilidades são equivalentes uma vez que resultam 
em probabilidades de erro de bit idênticas. 

Apresenta-se nas figuras 1 e 2 as codificações 
utilizadas no processo de simulação. Para índices de 
modulação superiores aos apresentados a constelação é 
construída de forma análoga. 


M=16 
M=4 

ua 3 Ata 1 oa | O ll 1 1 dd 1 
Di 11 1001 001 0101 1101 

s a +1 . * & . 
DO 10 1000 0000 0100 1100 

” . = | . Õ . “ 
- + | si 161 [8 au O O 1] O 1 1 Ú 
= 45 == | +1 + 3 

figura 1 Constelação e codificação Gray para 4-0AM e 16-0AM. 


M= 64 


A 111110 101110 100110 000110 001110 011110 010110 110110 

+7 . . “ . . e o . 
' 111010 101010 100010 000010 001010 011010 010010 110010 

+ s s . s . s s s 
; 111011 101011 100011 000011 001011 011011 010011 110011 

+53 » . Ê e e = - » 
111001 101001 100001 000001 001001 011001 010001 110001 

41 . e . . “ . . . 
111000 101000 100000 000000 001000 011000 010000 110000 

mi) & . s & & o Õ 
3 111100 101100 100100 000100 001100 011100 010100 110100 

a 1] a . ] [] ] ] E] 
- 111101 101101 100101 000101 001101 011101 010101 110101 

= "4 1] ] ] ] 1] LL] a O] 
” 111111 101111 100111 000111 001111 011111 010111 110111 

— | ) ] ] E] . = ] ] 
—7 — & —3 — +1 + + +7 

figura? Constelação e codificação Gray para 64-QAM. 
a 2 = ÇA * - 
5.2 Interferência intersimbólica. 


Probabilidade de erro de bit no 
receptor. 


Na presença de desvanecimento selectivo 
(únicamente) as componentes em fase e em quadratura à 
saida do desmodulador de M-QAM são dadas, 
respectivamente, por: |1| 
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u Ug)= 
Vo + + 49) 
oe po a Duda Ei gi ql 
uq tg)= 
V + o +00 (10) 


e L=VM; 

e Vorepresenta a amplitude do pico; 

º rei e rgi representa as amplitudes dos 
pré-cursores e pós-cursores.[1] 

No processo de simulação, depois de gerar uma 
sequência de bits de informação e os correspondentes 
símbolos em codificação Gray, determina-se para cada 
símbolo gerado, up(to) e ug(to), contabilizando o efeito da 
i.i.s. Sem perda de generalidade podemos fazer 
Vo(V2(L-1))=1 e definir: 


panda a (11) 
u tg)=u lt)! ro 


(12) 
u tg)=ultg ro 
Com esta normalização, no receptor têm-se os níveis t1, 
+3,..., +(L-1) possíveis para cada símbolo e os limites da 
zona de decisão são dados por 0, +t2, ...,t(L-2). 
Comparando up'(ty) e ug'(to) com a zona de decisão 
correspondente reconstitui-se o respectivo símbolo. 
Seguidamente, recorrendo à tabela de código de Gray, 
obtém-se os bits na recepção. 

Note-se que a modulação M-QAM exige um 
armazenamento prévio de u= logoM bits para obter um 
sinbolo. Assim, o número de bits a simular terá de ser 
divisível por u para se obter um número inteiro de 
simbolos. 


6 Determinação de uma assinatura por 
simulação computacional. 


O processo de determinação de uma Assinatura num 
sistema M-QAM, é constítuido pelos seguintes 
procedimentos: 

e Geração de bits; 

e Codificação em código de Gray; 

e Contabilização do efeito da i.i.s.; 

e Determinação da probabilidade de erro. 

A geração de bits permite obter os vários símbolos. 
Admite-se que, para efeitos de contabilização de 
interferência intersimbólica, um número de 4096 
símbolos é suficiente para definir uma probabilidade de 
erro de 10-3, sendo necessário gerar 4096-l0g9M bits para 
o efeito. 

Por leitura sequencial de u bits obtém-se os 
correspondentes símbolos com codificação de Gray. 
Seguidamente, contabiliza-se a i.i.s. em cada simbolo na 
forma: 


a4 


(13) 


N N (14) 
4, > L Ei Ei 2 ! Coil "ai 
— N — N 

onde N é o número de símbolos significativos para o 
cálculo de £, e &. Para efeitos de cálculo numérico, 
N=100 símbolos é suficiente para um cálculo utilizando 
reais com dupla precisão (64 bits). Esta "janela" implica 
a obtenção de 2:N símbolos adicionais de forma a que 
desde o primeiro ao último símbolo simulado esteja 
sempre presente i.i.s. na mesma forma, ou seja, l.i.s. quer 
devido a pós-cursores quer devido a pré-cursores. Este 
procedimento implica a geração de 2-N.logoM bits 
adicionais [5]. 

Uma vez contabilizada a i.i.s. em cada símbolo 
toma-se uma decisão baseada nos níveis de decisão 
indicados em $5.2 obtendo-se os símbolos no receptor. 
Seguidamente, por comparação com a tabela de 
codificação de Gray (Cf. fig's. 1 e 2), obtém-se os bits na 
recepção que permitem avaliar a probabilidade de erro. 

O passo seguinte consiste em repetir o processo atrás 
descrito formando o procedimento por lotes (Cf. 84). A 
escolha que se revelou mais adequada no procedimento 
por lotes foi n=10 e [= 10 que permite obter, na maioria 
das situações, amostras de P,p com correlação pouco 
significativa. À significância escolhida para o intervalo 
de confiança foi de 90% (a = 0.1) o que corresponde a t, 1. 
a/2= 1.383 na estatística t-Student para 9 graus de 
liberdade e uma significância de 90%. Quanto à 
correlação, considerou-se que as amostras eram 
incorrelacionadas quando |D| <0.3. 

O método descrito permite obter a probabilidade de 
erro de bit no receptor para um dado valor de b, fixando 
os parâmetros M, fp, à, 1, fo. Seguidamente, utilizando o 
método de Newton-Rapshon determina-se o valor de b 
que permite obter P,,=10-3 no receptor. O valor de b é 
determinado com uma precisão de 104. Finamente, 
variando fg é possivel obter b(fg) para P,p= 10-3, ou seja a 
Assinatura. 


7 Resultados obtidos. 


A limitação principal no processo descrito em 86 é o 
tempo de computação necessário à obtenção de uma 
Assinatura. No entanto, o principal objectivo deste 
processo é avaliar a precisão do método analítico descrito 
em [1]. Restará então comparar Assinaturas para um 
conjunto de parâmetros típico para avaliar a precisão 
correspondente. 

Devido ao elevado tempo de processamento não foram 
avaliadas Assinaturas para P,,= 10-89. Também foi 
excluido o caso M = 1024 por razões idênticas. 

Apresenta-se em seguida os resultados obtidos para 
M=4, 16, 64, 256 para o seguinte conjunto de 
parâmetros, típico num feixe hertziano digital de alta 
capacidade: f,= 140Mbit/s, a=0.5, 1=6.3 ns e para uma 
probabilidade de erro de 10-3 [5]. 


Técnica - 2/91 


7 Assinaturas em sistemas sem 


igualação 


Na tabela apresentam-se os valores numéricos e nas 
figuras 3 e 4 0 traçado das respectivas Assinaturas para 
um sistema sem igualação, utilizando o método analítico 
descrito em [1], assim como os valores da simulação por 
computador. 


Valores numéricos de Assinaturas sem igualação, para 


tabela | 
vários valores de M. 
5 — a 
+ ' Ê . é o 1 
a ; [ sã A dg 
ã | 
ss -— E | | 
E | | 
E dê | | 
= | | | 
| 
-35 | | 
ai | 
— qo | + ——>— we r A | us | 
-=40 — a dh [1] ED qi] 
fo [MHgzl 
=16 U=4 
figura3 Curvas Analíticas e Pontos de Simulação paraM=4e 
M= 16. 
B — = 
j p= UP ea, 
-10 || | T 
| 
-|5 || 
E =D + 
E du 
ah 
=) 
= 0) -— 0 no 1ó do sa 
fa |MHa! 
M=n4 M=7ba 
figura 4 Curvas Analíticas e Pontos de Simulação para M = 64 e 
M = 256. 
7.2 Assinaturas em sistemas com 


igualação LA E 5 


Utilizando um igualador simétrico com dois cursores 
(um pré-cursor e um pós-cursor) [1], obtiveram-se os 
valores numéricos das Assinaturas apresentadas na 
tabela 2 para P,,= 10-3 e 108, com igualação LAE. Nas 
figuras 5 e 6 apresentam-se os traçados gráficos. Os 
restantes parâmetros são apresentados em 87. 


5 DoInglês"Linear Adaptative Equalization” 
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CE OO RT 


tabela2 Valores Numéricos de Assinaturas com Igualação LAE, para 
M=ibeM=64. 
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figuraS Curvas Analíticas para P,, = 103 e 10% com Igualação LAE, 
para M= 16 e pontos de simulação. 
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figura 6 Curvas Analíticas para P,, = 107 e 10% com Igualação LAE, 
para M = 64 e pontos de simulação. 


Teo 

No processo analítico descrito em [1], existem duas 
aproximações inerentes ao método. À primeira diz 
respeito à majoração da 1.i.s. À segunda diz respeito ao 
número de símbolos interferentes para o fecho do olho, 
para uma dada probabilidade de erro, que se aproxima 
por Ni;zlogyP.p. Assim, no processo analítico, o 
problema reside em determinar o número inteiro 
mínimo de símbolos interferentes que permite obter um 
dado valor da probabilidade de erro que se aproxima pela 
expressão anterior. 

Para P,,=10-? apresenta-se de seguida a 
aproximação cometida no cálculo do número de símbolos 
interferentes, para vários índices de modulação. 

Como podemos verificar no quadro anterior, os valores 
de log; Pe, para M=4e 16 são praticamente coincidentes 
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com o número inteiro N; imediatamente superior. Para 
M=64 e 256 o inteiro imediatamente superior 
encontra-se mais afastado, sendo a diferença mais 
acentuada para M=64. Este efeito reproduz-se nos 
resultados obtidos por simulação (Cf. tab. 1, fig's. 3 e 4). 
Para M=4 e 16 os valores obtidos por simulação são 
praticamente coincidentes com os valores analíticos, 
sendo a majoração muito reduzida uma vez que a 
diferença anterior também é muito reduzida. Quanto aos 
valores para M= 64 e 256, a aproximação torna-se mais 
grosseira, como era de esperar, e a majoração da 
Assinatura é mais acentuada. Este efeito é mais notado 
para M= 64. 

As diferenças relativas da área sob as Assinaturas, 
expressas em unidades logarítmicas: 


= Analítico 
helio log 10 Simulação 


Mostra-se que o processo analítico constitui uma boa 
aproximação do processo simulado. 

Para a igualação LAÉE, com M= 16, verifica-se que os 
10 símbolos interferentes não são suficientes para 
majorar a Assinatura. Para esse facto contribui a 
pequena diferença existente entre o valor de logy P,, e 
N: conjuntamente com o facto da expressão N, = log; P.s 
ser aproximada. Na figura 3 podemos verificar que para 
N,/=20 (P,,=10-6) a curva se encontra bastante 
afastada demonstrando a influência de N; no andamento 
da Assinatura. Para M=64, a majoração é mais 
acentuada, verificando-se que os pontos de simulação 
tendem a aproximar-se mais da curva analítica, 
confirmando o referido para o caso de M=64 sem 
igualação. 

As diferenças relativas das áreas sob a Assinatura 
obtidas pelos dois processos são: 


Como conclusão, podemos afirmar que o processo de 
simulação confirmou com uma boa aproximação os 


resultados obtidos analíticamente, demonstrando que as 
diferenças entre os processos são muito reduzidas. 


8 Conclusão 


Desenvolveu-se um algoritmo de cálculo de 
Assinaturas de feixes hertzianos digitais na presença de 
desvanecimento selectivo. O algoritmo contempla a 
introdução da igualação LAE, desprezando-se possíveis 
imperfeições na implementação do sistema, obtendo-se 
Assinaturas para um caso ideal. O método possibilita a 
determinação de Assinaturas para os valores de 
probabilidade de erro de interesse no cálculo das 
margens definidas pela CCIR. 
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RESUMO 


O presente trabalho descreve um estudo numérico 
da análise de estabilidade da camada limite compres- 
sível e tridimensional sobre uma asa em flecha de 
avião, com vista à previsão do primeiro estágio da 
transição laminar-turbulenta, com especial relevância 
para tecnologias de escoamento laminar a usar nas 
próximas gerações de aviões. O método numérico pro- 
posto por Cebeci e Kaups para o cálculo de camadas li- 
mites laminares, compressíveis, sobre asas em flecha, 
é usado e testado. A estabilidade da camada limite, 
tridimensional, compressível, sobre uma asa em flecha 
é posteriormente investigada mediante uma análise de 
estabilidade linear e do método eN. Duas estratégias de 
integração do factor N são investigadas: para máxima 
taxa de amplificação, e para frequência e comprimento 
de onda constantes. Perturbações estacionárias e não-- 
-estacionárias são investigadas. Neste caso, os resulta- 
dos mostram que a amplificação de perturbações do 
tipo corrente cruzada é dominante em relação a insta- 
bilidades de Tollmien-Schlichting. 


1 Introdução 


É reconhecidamente aceite o facto de que as 
tecnologias do escoamento laminar apresentam um 
importante potencial, tendo em vista uma significativa 
redução na resistência aerodinâmica em cruzeiro de um 
avião. Manter o escoamento laminar sobre grande parte 
das asas e estabilizadores poderá conduzir a reduções na 
resistência aerodinâmica total da ordem dos 20%. 
Contudo, é igualmente aceite que este tipo de técnicas 
constituem uma tecnologia de alto risco, exigindo um 
empenhamento substancial na investigação a longo 
prazo. 

Técnicas modernas de produção permitem a obtenção 
de níveis de qualidade, em termos de acabamento de 
superfícies, compatíveis com os requisitos do projecto 
para escoamento laminar. Medições experimentais, por 
sua vez, revelaram que, sob certas condições, uma 
camada limite laminar pode subsisistir até 55% da 
corda, no caso de alguns aviões para utilizações 
vulgarizadas. Tais casos implicam, todavia, ângulos de 
flecha relativamente baixos e pequenos números de 
Reynolds (baseados na corda). Os modernos aviões de 
transporte, por sua vez, apresentam típicamente 
elevados números de Reynolds [0(107)-0(108)] bem 
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ABSTRACT 


The present work describes a numerical 
investigation based on linear stability analysis of the 
three-dimensional compressible boundary layer on a 
swept wing for the purpose of estimation of the first 
stage of laminar-turbulent transition. Appreciation is 
given to the concept of laminar flow technologies to be 
incorporated in the next generation of airframes. The 
numerical method proposed by Cebeci and Kaups is 
used and tested for solving the compressible flow 
around wings. The stability of the three-dimensional 
compressible boundary layer on a swept wing is later 
investigated using linear stability analysis in conjun- 
ction with the eN-method. Two different N-integration 
strategies are followed: maximum amplification rate, 
and fixed frequency and fixed wavelength of the 
disturbance. Travelling and stationary perturbations 
are investigated. In the present test case the results 
have shown that the amplification of crossflow-type 
disturbances prevail rather than Tollmien-Schlichting 
instabilities. 


como ângulos de flecha (do bordo de ataque) moderados 

(até 30º) combinados com escoamento transsónico. À 

tecnologia do escoamento laminar pode ser dividida em 

três ramos distintos: 

i) Recorre-se a um gradiente de pressão favorável para 
estabilizar a camada limite laminar (método 
passivo). Esta técnica ("Natural Laminar Flow") é 
aplicável desde pequenas aeronaves até aviões de 
médio porte; 

ii) Uma certa quantidade de sucção da camada limite 
laminar é utilizada, em conjunto com um gradiente 
de pressão favorável, para a estabilização da mesma 
(método activo). Esta técnica ("Hybrid Laminar Flow 
Control") permite a aplicação da tecnologia do 
escoamento laminar a aviões de transporte de médio 
porte com ângulos de flecha (do bordo de ataque) 
moderados; 

ili) À sucção da camada limite é utilizada para, por si só, 
permitir a estabilização da camada limite laminar 
(método activo). Esta técnica (LFC-"Laminar Flow 
Control") é indispensável para a manutenção de 
grandes regiões de escoamento laminar em aviões de 
transporte de grande porte, apresentando ângulos de 
flecha (do bordo de ataque) elevados. 

Assim, é inegável que a aplicação de técnicas de 
escoamento laminar a este tipo de aviões exigirá o desen- 
volvimento de métodos fiáveis de controlo (activo) da ca- 
mada limite compressível tridimensional. 
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No projecto de um perfil alar laminar para uma asa 
em flecha, é necessário considerar os diferentes mecanis- 
mos de instabilidade da camada limite: i)Instabilidades 
Tollmien-Schlichting (TS); ii)Instabilidades de corrente 
cruzada (CF); iii) Contaminação do escoamento laminar 
na zona de elevada curvatura do perfil (AL); iv) Vórtices 
de Górtler (GV). 

No presente estudo apenas são considerados os dois 
primeiros mecanismos, enquanto que aos restantes 
apenas se fará uma breve referência. 

A ideia de que a origem da turbulência proviria da 
instabilidade do escoamento laminar remonta ao séc. 
XIX. Uma revisão dos conhecimentos sobre instabi- 
lidades TS foi efectuada por Schlichting [1], contendo 
referências ao fenómeno de GV. Investigações remon- 
tando ao séc. XIX acerca da estabilidade da camada 
limite, foram confirmadas através das experiências clá- 
ssicas de Schubauer e Skramstad [2], dando crédito às 
teorias matemáticas desenvolvidas, para a situação 
bidimensional, por Lin [3]. O caso bidimensional foi 
estudado com grande detalhe nos anos 60 e 70 para os ca- 
sos particulares do escoamento de Blasius [4-10], e dos 
escoamentos Couette e Poiseuille [11], essencialmente 
por experimentação numérica. O desenvolvimento de 
novos métodos numéricos [12-14] veio permitir uma 
crescente eficiência na resolução numérica de equações 
tipo Orr-Sommerfeld. Estudos igualmente relevantes 
foram efectuados para o caso de escoamento não paralelo 
[6, 9], bem como a aplicação a escoamentos compres- 
síveis por Lees e Reshotko [15]. 

Para o caso de uma asa em flecha o escoamento é 
tridimensional, e este tipo de instabilidade viscosa é 
modificado devido à presença de escoamento de corrente 
cruzada, provocado pela flecha da asa. As instabilidades 
CF, dominantes em asas em flecha, são essencialmente 
semelhantes áquelas para o caso de um disco rotativo. 
Pelo facto de este último ser mais simples de estudar, 
foram efectuados trabalhos neste domínio (Malik et al. 
16] e Gregory et al., [17] ), considerando a influência de 
fôrças de Coriolis e da curvatura das linhas de corrente. 
Os efeitos de curvatura, tais como a curvatura das linhas 
de corrente e das superfícies, têm um efeito estabilizante 
no escoamento sobre superfícies convexas [18, 19]. O 
caso da estabilidade de escoamentos tridimensionais, 
incluindo instabilidades tridimensionais, tem sido 
igualmente estudado com detalhe[20, 21], bem como o 
efeito de sucção, arrefecimento, irregularidades 
superficiais e número de Mach, aplicado a asas em flecha 
ou a outros casos [22-25]. Importantes revisões de todos 
estes efeitos na estabilidade da camada limite foram 
apresentados por Mack [26] e Reshotko [271]. 
Instabilidaden da camada limite sobre a junção da asa à 
fuselagem, cu irregularidades no bordo de ataque, 
podem introduzir pequenos turbilhões, na camada limite 
laminar, que se propagam ao longo do bordo de ataque de 
asas em flecha. A transição por AL depende assim do 
ângulo de flecha, do raio do bordo de ataque, do número 
de Reynolds e da aceleração do escoamento na zona de 
elevada curvatura da superfície da asa. Um trabalho 
experimental de investigação deste fenómeno foi 
recentemente realizado por Arnal e Juillen [28]. No 
intradorso de perfis exibindo uma curvatura apreciável, 
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podem desenvolver-se vórtices de Górtler como 
consequência de instabilidades devidas a fôrças 
centrífugas em regiões de superfícies côncavas [29]. 
Quando o número de Reynolds excede determinado valor 
formam-se, nestes casos, na camada limite, pares de 
vórtices rodando em sentidos opostos. 

Investigações mais recentes têm focado as 
instabilidades secundárias na camada limite [30], e 
efeitos não lineares que caracterizam os últimos estágios 
do processo de transição [31-33], em virtude do aumento 
da capacidade de cálculo dos super-computadores. 
Contudo, no que se refere a códigos para previsão da 
transição, aqueles que presentemente têm maior 
utilização baseiam-se no método eN, um método 
semi-empírico que envolve o cálculo de taxas de 
amplificação de perturbações na camada limite, que são 
integradas para obter o chamado factor N. Este método 
assenta na premissa de que a localização da transição 
seja dada por um valor universal do factor N. Em 
camadas limites bidimensionais, onde as instabilidades 
TS são dominantes, o valor do factor N é 
aproximadamente 9 (e?=8000), para o local da 
transição. Tal não se verifica, em geral, para 
escoamentos tridimensionais, o que justifica ulteriores 
investigações. Contribuições para a análise numérica do 
problema [34, 35] foram importantes para o 
desenvolvimento de métodos bastante eficientes [36, 37], 
como SALLY (incompressível) ou COSAL 
(compressível). Estes dois métodos foram comparados 
entre si, por Malik e Orszag [38], bem como algumas das 
diferentes estratégias de integração possíveis em 
escoamentos tridimensionais, e com resultados de 
medições experimentais para asas incorporando LFC 
[39]. Os cálculos de estabilidade para o presente trabalho 
são efectuados com base num outro código, o COSALX 
[40, 41], uma versão modificada e melhorada do COSAL. 

Na indústria aeronáutica, o projecto de uma asa de um 
avião de transporte de longo curso reveste-se de 
particular importância, pelos objectivos descritos 
anteriormente. Um grande número de perfis alares é 
desenhado e vários testes preliminares são realizados 
(experimentação computacional). Destes, dois ou três 
perfis são seleccionados para experimentação em túnel 
aerodinâmico. À optimização de um destes conduzirá à 
forma definitiva. No caso de se pretender aplicar 
técnicas de LFC, a fase de escolha do perfil assume ainda 
uma maior relevância. É, portanto, indispensável dispór 
de um método relativamente barato do ponto de vista 
computacional, mas permitindo, simultaneamente, um 
cálculo suficientemente preciso da camada limite 
laminar. Um método com crescente utilização neste 
campo foi desenvolvido por Cebeci e Kaups [42]. Este 
permite efectuar o cálculo de camadas limites laminares, 
compressíveis sobre asas finitas em flecha com sucção. 
Por esta razão o método de Cebeci-Kaups foi escolhido 
para este trabalho. Contudo, algumas modificações 
foram incorporadas, tendo em vista uma mais 
apropriada geração dos perfis de camada limite, de modo 
a facilitar uma posterior análise de estabilidade baseada 
no método eN. A adequabilidade, para uma análise de es- 
tabilidade, do método Cebeci-Kaups, será investigada. 
Subsequentemente, proceder-se-á à análise de estabili- 
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dade da camada limite sobre uma asa em flecha, com 
base no método eN. Serão utilizadas as estratégias de 
integração para máxima taxa de amplificação da 
perturbação, e para perturbações com frequência e 
comprimento de onda constantes. As duas técnicas serão 
detalhadamente comparadas, e as suas vantagens 
discutidas. A contribuição de perturbações estacionárias 
para a transição da camada limite sobre uma asa em 
flecha será igualmente investigada. 


2 Modelo físico-matemático 


2.1 Camada limite laminar compressível 


sobre asas em flecha 


As equações que descrevem um escoamento de 
camada limite, compressível, são as equações de 
continuidade, as equações de Navier-Stokes e a equação 
da energia, assumindo a aproximação de camada de 
corte delgada. Quando escritas na forma estacionária, 
num referencial (z, 0, y), vide Figura 1, vêm: 


figura 1 Sistema de coordenadas e Geometria da asa planificada 


(x/c),, 


figura? Definição da Geometria do Perfil Alar 
e continuidade 


d d d ? 
— (, — — ( = (1) 
E (pzW) + o PU uid 0 
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sendo H a entalpia total e assumindo a aproximação 
9P/92=0. O ar é tratado como um gás perfeito, a 
viscosidade é obtida através da lei de Sutherland e o 
número de Prandtl é considerado constante. As 
condições de fronteira são dadas por: 


aH 
$=0: W=d Voy U=% (=) -0 (5a) 
u a 
Lil 
y>6: W=— W U->U H >»H (5b) 


Para ultrapassar o problema da singularidade ao 
longo da linha de estagnação, a equação (3) é derivada 
em ordem a 0, obtendo-se a partir de (1)(4) um novo 
conjunto de equações, impondo U =), aplicável ao longo 
desta linha [42]. 

O cálculo da camada limite é efectuado sobre a 
superfície planificada da asa, definida em função dos 
ângulos de flecha Ajg e Ark, para os bordos de ataque e 
fuga respectivamente. O escoamento de camada limite é 
analisado recorrendo a uma transformação de 
semelhança. Deste modo, a complexidade do problema é 
apreciávelmente reduzida, possibilitando que as 
equações que regem o escoamento tridimensional sejam 
escritas de um modo semelhante às de um escoamento 
bidimensional. Assume-se, fundamentalmente, que a 
asa planificada apresenta uma forma trapezoidal e que o 
gradiente de pressão ao longo da envergadura é 
desprezável. Assim, é possível eliminar a velocidade 
normal à superfície da asa como condição fronteira, a 
qual acopulava todas as variáveis, tornando-a um 
parâmetro das novas equações. Todavia, este tipo de 
análise não é aplicável às secções extremas da asa (raiz e 
bordo marginal). A semelhança da forma das secções, 
ignorando-se a torção da asa, bem como da distribuição 
de sucção (se considerada), são restrições impostas pelo 
presente método. 

A introdução da variável da semelhança n e de um 
potencial vector (!P, D), definidos por: 


W 12 
do=( e aus (6) 
 P,H,2 
12 
o W 
v=2(p 0 W) f(n,0) , f=— (Va 
é W 
Ê 
12 U 
o U (7) 
0=22[puW,) vig(n,0) =— « 
e Em E w, h ) ' 5 U, 


permite tornar as equações que regem o escoamento 
independentes de z, podendo a solução ser obtida na 
direcção 6, a partir do bordo de ataque. 
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As equações resultantes da transformação poderão 
ser encontradas em |42|. 


2.1.1 Modelo numérico 

Para a resolução das equações de camada limite 
resultantes é utilizado, na sua discretização, o método de 
diferenças finitas conhecido por "Keller box method”, o 
qual é adequado à resolução de equações às derivadas 
parciais parabólicas, obtendo-se uma precisão de O(h2). 
Detalhes deste método poderão ser encontrados em [43, 
44]. 

A discretização é efectuada numa malha (vide Figura 
3), cujo espaçamento entre nós é dado por uma 
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figura3 Malha utilizada no cálculo da camada limite. Factor de 
ampliação da dimensão perpendicular à superficie: 100. 

progressão geométrica de razão constante, a qual é 

adequada para o fim último de análise de estabilidade. 


23 Estabilidade linear da camada limite 


tridimensional compressível 


A derivação das equações para a análise de estabili- 
dade linear [26] é efectuada a partir das equações de 
continuidade, Navier-Stokes, e energia para um fluido 
viscoso, condutor de calor, tomado como um gás perfeito, 
num sistema de coordenadas cartesianas. 

As quantidades instantâneas são decompostas nas 
suas componentes referentes ao campo médio e 
flutuações não-estacionárias: U=U +u, e análogamente 
para V, W, T,P,p, wu, À, e k Assumindo escoamento 
quase-paralelo: 


U=U) : T=T0) 
V=0 : P=P (X,2) (8) 
W=W6y) 


as variáveis instantâneas são substituidas nas equações 
de continuidade, quantidade de movimento, energia e 
estado pela sua forma decomposta e, seguidamente, os 
termos apenas referentes ao campo médio são subtraídos 
às equações resultantes, originando as equações para as 
flutuações, envolvendo u, v, w, O, p. Os coeficientes de 


viscosidade e condutividade térmica são considerados 
apenas função da temperatura, pelo que as suas 
correspondentes flutuações podem ser escritas como: 


dA dk (9) 


As equações resultantes são adimensionalizadas em 
relação a U, e 6" (definição compressivel ou 
incompressivel), e aos valores exteriores à camada limite 
das variáveis de estado; os coeficientes de viscosidade são 
tornados adimensionais através de nu, e a condutibilidade 
térmica através de Cone. Às equações são linearizadas, 
ou seja, termos quadráticos nas flutuações são 
desprezados, vindo então: 

€ continuidade: 
+ pol + U- + wir = () 
dz 


dr (= du =) 
dy dx 


—+p|—+— +— 
dt dx dy dz 


(10) 


e quantidade de movimento segundo x: 


-=(du —=du aU mm () Lj 
p(—+ O W— )= 
dt dx dy dz 


2 Z É 
du du du d w 
ar 3y Em didy dxdz (1 1) 


i x E, E 

E (ta au du do w dudT “du du' 

A-p( = +—+ E (4) 
e dxdy aTrdy toy dx. 


E p= mi 
dU ad dudTdU. | 
+——— | + 
E A , à sd À p , 
dT a dy dy dTº dy dy 


€ quantidade de movimento segundo y: 


= duo— hi) — ph qi 
o( —+U—+W— |= 
Ot dx dz 


A 2 » E. dh 
-=tyu du du du do tu 
dr 


— + — + 
2 1. 2 2 
dy dx dz 


didy dydz 


(12) 
dy 

du dTov 3U ab dW al 

(2 +++ 


dy dx 


+ 


dT'* dyày dy dz: 


Ri Fuad du 3] 
34dT dT 


dy *dx dy dz 


e quantidade de movimento segundo z: 


dW Wo) 
2 2 2 à 
= w p(— du du =; 


—— + — + 
2º . dxdz 
dy dx 


dydz 


- cc —= [O Ww 5 e A 
ad p( + er (13) 
K 2 dGydz dxdz 

E (0 =) 
== a [o O e 
dTdy dy dz 


a To 4 dy 


dWob ) 
= |+ —— — 


De 
d HdTdW , 
ar ge Upa 
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e energia: 
a ar o” 90) 
(= + EE + W— )= 

dé dx dy dz 


va [30 2%0 9% idkd'T 
PrR ax? ay? az? kdT dy" 
(14) 


jm ce e 


adEdTo0 1 dk ea | 
=—— += (— 6 |+ 
kdTdydy k aT? dy 


ytr- MM 
+——+|3 
R 


su | 2 
agi E EA) 
dy Xdz dy dT* dy dT* dy. 


(o -) 
— + — |) + 


dy Xdôy dx 


€ equação de estado: 
rr. B 


p==+= 


(15) 


b 
Estas equações admitem soluções da forma 
harmónica: 


Lu, U, W,p,r, 9! = 


(16) 
[0(y), v(y), uMy), ply), rly), 1" exp| il Jadx + Jfdz-wt) | 
com condições de fronteira, por exemplo: 
y=0: u,v,w,r,8=0 
y>o; u(y), v(y), w(y), r(y), O6)=0 (17) 


As equações finais [26] resultam da substituição das 
expressões (16)-(17) nas equações (10)-(15). 

As funções Q(y), d(y), W(y), Pl(y), Ply) e O(y), são 
usualmente designadas por funções próprias ou funções 
de amplitude, na medida em que descrevem a 
distribuição da amplitude das perturbações através da 
camada limite. As soluções harmónicas consideradas 
apresentam, no caso geral, números de onda nas 
direcções x e z, respectivamente a e Ê, e frequência w, 
complexos. No entanto, caso sejam reais, a perturbação 
será estável, propagando-se no plano (x, z) com 
velocidade e amplitude máxima constante (estabilidade 
neutral). Pelo contrário, se a, 8, w forem complexos, a 
amplitude da perturbação variará à medida que esta se 
propaga. À procura deste tipo de soluções resulta num 
problema de valores próprios w =w(a, B) envolvendo seis 
quantidades. Contudo, apenas duas destas podem ser 
determinadas da solução do problema de valores 
próprios, pelo que as restantes terão de ser especificadas 
de alguma forma. Assim, surgem duas teorias: 

i) Teoria da Amplificação Temporal: 


u= O(y).explco, £). exp 


( fadx+ fBde-a,4) 
a=a , B=b. 


v=uw tim, 
r L Fr 
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módulo do vector número de onda 
O qa 
el=(02+p2) 
ângulo entre o vector número de onda e a direcção x 


na (2) 
a 


A amplitude da perturbação variará no tempo com uma 
taxa w;, tal que: 


ldA | 

Ade 
w;<O - aperturbação é amortecida 
w;=0 -— a perturbação é estável 
w;>0 -— aperturbação é amplificada 


ii) Teoria da Amplificação Espacial: 

u=0 (y.exp| —(Ja, dx +18, de ) |expli(Sa, dx+18, de= 1) | 
B=B +16, 

A perturbação será amplificada (ou amortecida) nas 
direcções x e/ou z dependendo dos valores de a; e B;. 


Apesar de esta teoria ter um maior significado físico do 
que a anterior, implica uma análise mais complexa. 


(= q) , q=a +ia. 
r r i 


2.2.1 Modelo numérico 

O método utilizado neste trabalho, o COSALX [40,41], 
permite a análise da estabilidade linear de camadas 
limites tridimensionais. Os perfis de velocidade e 
temperatura da camada limite poderão ser gerados pelo 
método proposto por Cebeci e Kaups [42]. 

Pela razão apontada em 2.2, é seguida a teoria da 
amplificação temporal. No entanto, é utilizada a 
transformação de Gaster [35], a qual permite efectuar a 
conversão da amplificação temporal na sua 
correspondente espacial: 


d d 
— =|Re(V )| — (18) 
rd dx, 

sendo Ve =(dw/da, dw/9B) o vector grupo-velocidade. 
Assim, resulta a amplificação espacial na direcção do 
vector grupo-velocidade: 


aee ps (19) 
— *g Re(V q) 


As equações (10)(17) podem ser escritas na forma 
matricial [36, 37): 
E (20) 
AD'p+BDP+C.P=0 
sendo discretizadas, numa malha desfasada 
perpendicular à superficie, num esquema de diferenças 
finitas com uma precisão de O(h2). 
O problema de valores próprios, introduzido em 2.2, é 
resolvido da seguinte forma: 
1) Não havendo uma boa aproximação inicial para o 
valor próprio, é resolvido o problema generalizado, ou 
seja: 


[B!A-w1)4=0 (21) 
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sendo w o valor próprio (complexo) e | a matriz 
identidade. É utilizado um método espectral, o qual é 
globalmente convergente; 
ii) Havendo uma boa aproximação inicial para o valor 
próprio, é utilizado um método de convergência 
cúbica, baseado num processo iterativo ("Inverse 
Rayleigh Iteration”), para melhorar aquele valor, ou 


seja: 


(22) 


[DA a q 


Qk+h) (k+ 1) 
|v 4 


A análise de estabilidade linear permite a utilização 
do método eN na previsão da transição do escoamento 
laminar para turbulento. Como se referiu 
anteriormente, este método consiste essencialmente na 
integração das taxas de amplificação da perturbação no 
escoamento. Em escoamentos tridimensionais a referida 
integração pode ser efectuada segundo diferentes 
estratégias. De entre as várias opções possíveis apenas 
se fará referência áquelas utilizadas no decorrer do 
presente trabalho: 


i) Máxima Taxa de Amplificação 
Impondo inicialmente a orientação da perturbação 
relativamente ao escoamento exterior, Wg, € O seu 
comprimento de onda, À,,, é procurada, para a 
estação inicial, uma perturbação amplificada 
(w;>0). Quando esta é encontrada, a integração é 
efectuada sucessivamente (até à estação final, ou até 
se encontrar uma região estável) para as 
perturbações, com a frequência pretendida, 
apresentando a maior taxa de amplificação em cada 
estação. Tal é conseguido através do 
grupo-velocidade; 

11) Frequência constante e Comprimento de Onda 

Constante 

Impondo inicialmente py e fixando À,, é procurada, 
para a estação inicial, uma perturbação amplificada. 
Quando esta é encontrada, são efectuadas iterações, 
alterando a orientação da perturbação, por forma a 
que a perturbação (instável) apresente não só o À, 
inicialmente imposto como também a frequência 
pretendida. Esta é, então, seguida na integração laté 
à estação final, ou até se encontrar uma região 
estável). 


As diferentes estratégias de integração do factor N 
advêm do facto de que este não se encontra 
completamente definido, sem que se imponha a forma da 
perturbação a seguir em cada ponto, sobre a asa, e se 
defina o caminho de integração. Foi referido em 2.2, que 
a relação que exprime o problema de valores próprios 
envolve três quantidades complexas a, 5 e «, fornecendo 
apenas duas relações reais entre estes parâmetros. Para. 
o caso em questão, a Teoria da Estabilidade Temporal, 
impõe-se, como igualmente se mencionou, que a e É 
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sejam reais, o que fornece duas condições adicionais, 
Restam então, dois parâmetros arbitrários entre a,, br, 
we w;. Esta arbitrariedade pode, assim, ser resolvida de 
várias maneiras, mas é igualmente importante que a 
solução não seja desprovida de significado físico. 

Então, o factor N é calculado da seguinte forma: 


A EM 
N=in( O )= | ES ds 
A! JsIRe(V)| 
o u g 


3 Resultados 


(23) 


Com o objectivo de testar as possibilidades do presente 
método foi efectuada, numa primeira fase, um cálculo 
para o extradorso de uma asa em flecha. Trata-se de uma 
asa infinita, apresentando uma flecha Ag =Arp=18º, 
num escoamento a número de Mach, Mw=0.7. 
Prescreveram-se as condições iniciais e fronteira 
necessárias à aplicação do método, nomeadamente a 
descrição da geometria (y/c)p us (x/c)y, à distribuição de 
Cp e parâmetros do escoamento não perturbado, tais 
como Ps, Te e Pr. O cálculo da camada limite permite 
obter as componentes radial e circunferencial da 
velocidade invíscida, bem como os perfis de 
velocidade/temperatura e suas primeira e segunda 
derivadas, necessários para uma posterior análise de es- 
tabilidade. Com estes valores são calculados os 
parâmetros integrais de camada limite, espessura de 
deslocamento e espessura de momento, assim como os 
coeficientes de tensão de corte superficial nas direcções 
circunferencial e radial. De todas estas grandezas, para a 
análise de estabilidade, ver equações (10)-(15), só 
interessa reter os perfis de velocidade e temperatura e 
suas derivadas, e a espessura de deslocamento como 
escala de adimensionalização. 

As Figuras 4a) e 4b) mostram, respectivamente, a 


iii a 
= 
[as 


figura 4a Curvas de nível de numero de Mach através da Camada 
Limite. Factor de Ampliação da Dimensão Perpendicular à 
Superficie:100 
distribuição de curvas de nivel de número de Mach e 
temperatura, sobre o extradorso da asa, onde é notório o 
principio de separação da camada limite para valores 
apreciáveis da corda. As figuras 5a) e 5b) reportam os 
coeficientes de tensão de corte superficial 
respectivamente nas direcções circunferencial e radial 
sobre o extradorso da asa. Note-se que a separação se 
inicia para os perfis de velocidade na direcção 
circunferencial, o que está de acordo com as Figuras 4a) e 
4b).A título de exemplo, apresentam-se também perfis 
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figura 4b Curvas de Nível de Temperatura através da Camada Limite. 
Temperatura em Kelvin. Factor de Ampliação da Dimensão 
Perpendicular à superficie: 100. 
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figura 5a Evolução do Coeficiente de Tensão de Corte Superficial, na 
Direcção Circunferencial, ao longo da corda 
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figura Sb Evolução do Coeficiente de Tensão de Corte Superficial, na 
Direcção Radial, ao longo da corda 


adimensionais de velocidade circunferencial e radial e de 
temperatura e suas derivadas (vide Figuras 6a)-6c)), 
para diferentes estações sobre o extradorso da asa, 
utilizados na subsequente análise de estabilidade. 

Na medida em que se trata de uma asa em flecha, 
apresentando regiões de elevado gradiente de pressão ao 
longo da corda (favorável junto ao bordo de ataque), 
verifica-se o desenvolvimento de um perfil de velocidade 
de corrente cruzada. Na região do bordo de ataque, as 
linhas de corrente são deflectidas devido a gradientes de 
pressão na direcção da envergadura. Na camada limite o 
fluido é retardado junto à superfície, sendo assim ainda 
mais deflectido, dando origem a um perfil de velocidade 
perpendicular à direcção das linhas de corrente. Este 
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figura 6a Perfil de velocidade Circunferencial e Derivadas, a 
X/€ = 0.00004 sobre o Extradorso 
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figura 6c Perfil de Temperaturas e Derivadas, a X/C = 0.2995 sobre o 
Extradorso 


tipo de perfil (vide Figura 7) apresenta ponto de inflexão, 
pelo que é dinâmicamente instável. De facto, apenas 
uma pequena quantidade de escoamento de corrente 
cruzada (“crossflow”, CF) é necessária para que se 
verifiquem instabilidades deste tipo. 

Além de instabilidades CF, as quais são 
características de camadas limites tridimensionais, 
verificam-se igualmente instabilidades do tipo invíscido 
(relacionadas com perfis exibindo ponto de inflexão) e 
viscoso mais semelhantes áquelas que surgem em 
camadas limites bidimensionais, embora modificadas 
devido à presença de escoamento CF. Este tipo de 
instabilidade, usualmente designada por 
Tollmien-Schlichting (TS) ou "streamwise”", engloba 
todas as instabilidades que não estejam associadas a 
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DIRECÇÃO DO ESCOAMENTO 
POTENCIAL 


Y wW 


U PONTO DE INFLEXÃO 


CORRENTE CRUZADA 


figura7 Representação esquemática dos perfis de velocidade 


escoamento CF ou à inversão do escoamento. 
Típicamente, para uma asa em flecha, a região de 
instabilidades do tipo TS surge mais tarde do que a das 
do tipo CF. 

A Figura 8 mostra as previsões das taxas de 
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figura8B Taxas de Amplificação Temporal em função da Frequência 

da Perturbação, para diferentes Comprimentos de Onda, a 

X/C = 0.1223 sobre o Intradorso 
amplificação (temporais) verificadas para diversos 
comprimentos de onda de perturbações, em função da 
frequência destas. Refere-se a X/C=0.1223, sobre o 
intradorso do perfil, o que permite observar 
perturbações quer do tipo CF, quer do tipo TS. As 
instabilidades do tipo CF são claramente dominantes, 
verificando-se para menores comprimentos de onda, mas 
não tão reduzidos como para as primeiras estações, € 
atingindo as máximas taxas de amplificação mais ou 
menos próximo de pg= 88º. Estas são obtidas para ondas 
com uma frequência de aproximadamente 1 kHz, 
embora as ondas de frequência zero (perturbações 
estacionárias) apresentem também elevada 
amplificação, correspondendo mesmo à máxima para 
Aw=0.001. Felevante é igualmente a curva para 
Aw= 0.004, apresentando dois máximos, e sendo já, para 
este comprimento de onda, mais importantes as 
instabilidades do tipo TS. Note-se também que se 
verifica amplificação das perturbações, ainda para o 
mesmo comprimento de onda, com uma orientação em 
relação ao escoamento exterior pyg<O0, mas para 
menores frequências do que para yg=0. Para esta 
perturbação a máxima taxa de amplificação, na referida 
estação, é atingida para wy;=38º, uma perturbação do 


tipo TS oblíqua. Uma importante conclusão, que se 
retira também da Figura 8, é a de que a banda de 
frequências possíveis para uma perturbação alinhada 
com o escoamento exterior é mais alargada do que aquela 
para uma perturbação com uma orientação próxima da 
direcção típica CF. Estes resultados, não só para esta 
estação específica, mas ao longo de toda a corda, 
constituem uma indicação de grande relevância para a 
pesquisa de perturbações relevantes no processo de 
transição, atendendo a que não se possuem à priori 
estimadores das suas frequências e comprimentos de 
onda, competindo ao investigador descobri-las através de 
um cálculo interactivo. 

A investigação da transição é efectuada, como já se 
referiu, recorrendo a um método semi-empírico, o método 
eN, em conjunto com a teoria linear da estabilidade. De 
acordo com o mencionado em 2.2.1, uma das hipóteses, 
para proceder à integração para a obtenção do factor N, 
consiste em integrar sucessivamente as taxas de 
amplificação máxima para uma dada frequência. A 
Figura 9 mostra os factores N, para máxima taxa de am- 
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figura9 Evolução do Factor N, para perturbações de diversas 


frequências, pelo Método da Máxima Taxa de 
Amplificação, sobre o Extradorso. 


plificação, abrangendo uma larga banda de frequências 
das perturbações, para o extradorso do perfil. Resultados 
semelhantes foram obtidos para o intradorso. As 
integrações são interrompidas ao atingir-se a região 
estável. Note-se como as instabilidades com uma 
frequência mais elevada rápidamente se tornam estáveis 
à medida que nos afastamos da região de elevada 
curvatura. Os valores mostram-se elevados a pouca 
distância do bordo de ataque, De facto, admitindo válido 
o critério de que a ocorrência de transição se verifica 
para N=9, esta verificar-se-ia apenas a 4-6% da corda. 
Tal critério, contudo, não é universal. Além disso, o 
método de seguir uma perturbação periódica no tempo e 
não no espaço, constitui uma estratégia de integração do 
factor N desprovida de significado físico. 

Para o caso de perturbações de frequência é 
comprimento de onda constantes é significativo 
estabelecer a diferença entre aquelas que apresentam 
uma orientação típica de CF, e as restantes designadas 
por TS. Para estas últimas, as perturbações 
significativamente amplificadas encontradas são 
escassas, apresentando comprimentos de onda 
relativamente mais elevados, e iniciando a amplificação 
bastante tarde. Perturbações exibindo maiores 
comprimentos de onda iniciam a sua amplificação mais a 
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juzante. Estas evoluções estão representadas na Figura 
10, para o extradorso do perfil. A banda de perturbações 
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figura 10 Evolução do factor N, para perturbações do tipo 
Tolimien-Schlichting de Frequência e Comprimento de 
Onda Constantes, sobre o Extradorso. 
amplificadas do tipo CF é bastante mais interessante. Os 
resultados de uma investigação mais intensiva estão 
expressos na Figura 11 para o intradorso do perfil. 
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figura 11 Evolução do factor N, para perturbações do tipo Corrente 
Cruzada de Frequência e Comprimento de Onda 
constantes, sobre o Intradorso. 


Destaque-se que, para comprimentos de onda 
suficientemente baixos, ou seja 0.0003< Ay <0.0005, as 
taxas de amplificação são elevadas, aproximando-se 
mesmo dos valores obtidos através do método da máxima 
taxa de amplificação. A dependência da frequência não é 
particularmente relevante, parecendo apenas que uma 
perturbação com uma frequência mais elevada se torna 
estável mais cedo, a exemplo do verificado para o método 
da máxima taxa de amplificação, embora esse efeito não 
seja agora tão evidente. De facto, para perturbações 
dentro da gama de comprimentos de onda acima 
referida, a região estável é invariávelmente atingida 
muito perto do bordo de ataque, quedando-se o factor N 
por valores sempre inferiores a 5 para o intradorso e 8 
para o extradorso. Importante é a dependência, já 
observada para TS, em relação ao comprimento de onda. 
Perturbações de comprimento de onda superior aos 
acima referidos iniciam a amplificação claramente mais 
longe do bordo de ataque, atingem a região estável 
igualmente mais longe (para os maiores Ay investigados 
não chegam mesmo a atingi-la), mas apresentam taxas 
de amplificação muito inferiores, implicando valores do 
factor N relevantes apenas a 50 60% da corda. 

Assim, estes valores do factor N foram considerados 
situarem-se abaixo dos limites adequados. 


J.C.F. Pereira e J.M.M. Sousa 


4 Conclusões 


Uma análise baseada na teoria linear da estabilidade 
do escoamento na camada limite tridimensional 
compressivel sobre uma asa em flecha de um avião a 
Mach igual a 0.7 foi resumidamente apresentada e 
discutida. 

Duas estratégias diferenciadas de integração do factor 
N, máxima taxa de amplificação e frequência e 
comprimento de onda constantes, foram extensivamente 
comparadas. O método da máxima taxa de amplificação 
produz valores muito elevados do factor N ainda próximo 
da região do bordo de ataque da asa em flecha. A 
estratégia de integração para frequência e comprimento 
de onda constantes produz valores do factor N bastante 
mais reduzidos. Esta última apresenta, no entanto, um 
maior significado físico. Uma investigação da transição 
utilizando a referida estratégia de integração é mais 
complexa, uma vez que se verifica a possibilidade de 
seguir uma perturbação com uma orientação típica de 
corrente cruzada ou Tollmiem-Schlichting, bem como 
uma significativa dependência do valor do comprimento 
de onda. Contudo, para uma asa em flecha, as 
perturbações do tipo corrente cruzada são, em geral, 
dominantes. 
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